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Chapter  1 
INTRODUCTION 


Assuming  a  Gaussian  input  signal,  design  of  the  optimum  radar  processor  is 
well  known.  Specifically,  consider  a  coherent  train  of  N  received  pulses,  with  ran¬ 
dom  initial  phase  and  amplitude,  embedded  in  clutter  plus  noise.  Let  the  clutter 
plus  noise  be  referred  to  as  the  disturbance  process.  For  Gaussian  inputs,  it  is 
known  that  the  probability  of  detection  is  maximized  for  a  constant  false  alarm 
probability  by  maximizing  the  signal-to-disturbance  ratio  (SDR)  [1],  To  maximize 
SDR,  a  transversal  filter  is  placed  before  the  likelihood  ratio  test  (LRT)  detector. 
The  optimum  transversal  filter  is  completely  described  by  a  set  of  N  complex 
weights  whose  values  depend  upon  the  disturbance  statistics  and  the  target  sig¬ 
nal  model  [2]. 

Analyses  of  radar  performance  usually  proceed  by  assuming  a  linear  receiv¬ 
er  [3].  In  the  presence  of  strong  signals,  however,  actual  receivers  are  subject  to 
saturation  effects.  This  nonlinear  behavior  can  cause  a  significant  change  in  the 
statistics  of  the  disturbance.  The  goal  of  this  study  is  to  determine  the  amount  of 
degradation  in  the  signal-to-disturbance  ratio  at  the  transversal  filter  output  when 
a  nonlinearity,  whose  input-output  characteristic  models  the  receiver  saturation 
curve,  is  included  within  the  RF  amplifier. 

Three  different  situations  are  investigated  as  indicated  below: 

Case  I  -  The  radar  receiver  is  modeled  as  being  linear  and  the  weights  of  the 
transversal  filter  are  chosen  so  as  to  maximize  the  SDR. 

Case  II  -  The  radar  receiver  is  modeled  as  being  nonlinear.  Even  so,  the 
weights  of  the  transversal  filter  are  set  equal  to  those  used  in  case  I.  This  case 
can  arise  when  the  design  engineer  is  either  unaware  of  the  nonlinearity  or  elects 
not  to  compensate  for  it. 

Case  in  -  The  radar  receiver  is  modeled  as  being  nonlinear.  However,  now 
the  weights  of  the  transversal  filter  are  selected  to  maximize  the  SDR  taking  into 
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account  the  receiver  nonlinearity. 

The  amount  of  degradation  in  SDR  caused  by  ignoring  the  nonlinearity  is  ob¬ 
tained  by  comparing  cases  I  and  II.  Comparison  of  cases  n  and  III  reveals  the  im¬ 
provement  to  be  gained  by  accounting  for  the  nonlinearity  .  Finally,  the  loss  due  to 
the  nonlinearity  even  with  the  optimum  transversal  filter  is  evaluated  by  comparing 
cases  I  and  III. 

In  chapter  2  the  analytical  model  used  to  analyze  the  radar  receiver  is  pre¬ 
sented.  Analytical  expressions  are  derived  in  chapter  3  for  the  nonlinear  respons¬ 
es  at  each  point  of  the  receiver  model.  The  relations  between  autocorrelation 
functions  at  the  various  receiver  stages  are  presented  in  chapter  4.  The  autocor¬ 
relation  function  at  the  transversal  filter  input  is  derived  in  chapter  5  for  the  Gaus¬ 
sian  case  where  specific  models  are  assumed  for  the  target  signal,  clutter,  and 
noise.  Expressions  for  the  signal-to-disturbance  ratio  at  the  transversal  filter  out¬ 
put  are  obtained  in  chapter  6  for  cases  I,  n,  and  HI.  Up  to  this  point,  all  of  the  anal¬ 
ysis  is  carried  out  in  the  time  domain.  A  frequency  domain  interpretation  is  given 
in  chapter  7.  Finally,  computer  generated  results  are  presented  in  chapter  8. 
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Chapter  2 
RECEIVER  MODEL 


The  model  used  to  analyse  the  receiver  is  illustrated  in  Fig.  2.1 .  Only  the  sig¬ 
nal  processing  stages  before  the  LRT  detector  are  included  in  the  figure,  let  rk(t) 
denote  the  signal  at  a  point  k. 


Figure  2. 1  -  Receiver  Model. 

As  shown  in  Fig.  2.1 ,  the  received  signal  r^t)  is  first  amplified  and  filtered  by 
the  RF  amplifier.  The  output  r3(t)  is  then  multiplied  by  the  local  oscillator  signal 
and  low  pass  filtered  to  produce  the  baseband  signal  r5(t).  The  range  gate  selects 
the  successive  pulse  returns  from  a  particular  range  cell.  These,  in  turn,  are  pro¬ 
cessed  by  a  bank  of  transversal  filters  so  as  to  maximize  the  signal-to-disturbance 
ratio  at  the  input  of  the  LRT  detector.  Each  transversal  filter  is  tuned  to  a  different 
Doppler  frequency.  The  transversal  filter  whose  Doppler  frequency  coincides  with 
that  of  the  desired  signal  yields  the  maximum  output  in  the  absence  of  distur¬ 
bance  and  is  the  only  one  analyzed  in  this  work.  In  this  chapter  the  analytical 
model  for  each  block  is  discussed  in  detail. 
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For  simplicity  of  analysis,  the  receiver  nonlinear  behavior  is  assumed  to  be 
confined  to  the  RF  amplifier.  As  seen  in  Fig.  2.2,  the  RF  amplifier  is  modeled  as  a 
zero-memory  nonlinear  amplifier  cascaded  with  an  ideal  bandpass  filter  having 
unity  gain.  Observe  that  the  nonlinear  amplifier  output  is  denoted  by  r2(t). 


Fiaure  2.2  -  Block  Diaaram  of  the  RF  Amplifier. 


2.1.1  The  Nonlinear  Amplifier: 


Receiver  saturation  effects  are  accounted  for  by  the  zero-memory  nonlinear 
amplifier  whose  input-output  characteristic  is  given  by 


r2(t) 


(2.1) 


Thus,  r2(t)  is  propotional  to  the  integral  of  a  zero-mean  Gaussian  density  function 
with  variance  ag2.  A  sketch  of  r2(t)  versus  r-|(t)  is  shown  in  Fig.  2.3.  The  amplifier 
has  zero  memory  in  the  sense  that  the  present  value  of  the  output  is  determined 
only  by  the  present  value  of  the  input.  Observe  that  L  represents  the  saturation 
level  of  the  non-linearity. 
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The  slope  of  the  input-output  characteristic  is 


dr2(t) 

d?7(t) 


f 


v 


r?  (t) ' 

2o2 


g  j 


(2.2) 


r2(0 


Figure  2.3  -  Input-Output  Characteristic  of  Zero-  Memory  Nonlinear 

Amplifier. 

Let  the  linear  gain  of  the  RF  amplifier,  denoted  by  G,  be  defined  as  the  slope  of 
the  input-output  characteristic  evaluated  at  the  origin.  It  follows  that 


dr2(t) 
dr,  (t) 


r,  (t)  =  o 


2L 


(2.3) 


For  small  enough  inputs,  the  nonlinear  amplifier  output  can  be  determined  from 
the  linear  characteristic 


r2(t)  =  G  r,(t)  . 


(2.4) 
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Note  that  the  linear  gain  involves  the  ratio  of  L  to  og  .  When  r-, 

(2.1)  and  (2.4)  yield 

0.6826L  ;  Eq(2.1) 

0.7979L  ;  Eq(2.4)  • 

Hence,  when  r-|=ag,  the  ratio  (in  dB)  of  the  outputs  from  the  nonlinear  and  linear 
characteristic  is 

20L°9ioS>  -I-355*6  '  <2-6> 

Conventionally,  the  1  dB  compression  point  determines  the  linear  region  of  a  non¬ 
linear  amplifier.  It  follows  that  the  nonlinearity  of  Eq.(2.1)  can  be  considered  to  be 
approximately  linear  over  the  range 

(2-7) 


=ag,  Equations 

(2.5) 


2.1.2  The  Bandpass  Filter 

Let  fc  denote  the  carrier  frequency,  in  Hertz,  of  the  received  signal.  The 
magnitude  of  the  ideal  bandpass  filter  transfer  function  is  shown  in  Fig. 2. 4.  Note 


Wmm 


Figure  2.4  -  Transfer  Function  Magnitude  of  Ideal  Bandpass  Filter. 
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that  the  filter  has  unity  gain  over  the  entire  bandpass  region  and  zero  gain  other¬ 
wise.  The  filter  bandwidth  Brf  is  the  smallest  bandwidth  required  to  pass  undis¬ 
torted  the  desired  signal  portion  of  r2(t).  Note  that  the  filter  removes  the  nonlinear 
products  in  r2(t)  centered  at  harmonics  of  the  carrier  frequency  along  with  compo¬ 
nents  of  the  noise  which  fall  outside  of  the  desired  signal  bandwidth. 

2.2  The  Mixer 

The  mixer  consists  of  the  multiplier,  local  oscillator,  and  lowpass  filter.  It 
serves  to  translate  the  spectrum  of  r3(t)  from  ±fc  to  dc.  Assuming  the  local  oscilla¬ 
tor  signal  is  a  unit  amplitude  sinusoid  at  the  carrier  frequency  fc,  the  input  to  the 
lowpass  filter  is 


r4  (t)  =  r3(t)cos(27dct)  .  (2.8) 

The  lowpass  filter  is  assumed  to  be  ideal  with  a  transfer  function  magnitude  as 
shown  in  Fig.  2.5.  Once  again,  the  filter  gain  is  unity  over  the  passband  and  zero 
otherwise.The  filter  bandwidth  BL  is  the  smallest  bandwidth  required  to  pass  un¬ 
distorted  the  desired  signal  portion  of  r4(t).  Note  that  r5(t)  is  a  baseband  signal. 

I  HL(f) ! 


Figure  2.5  -  Transfer  Function  Magnitude  of  Ideal  Lowpass  Filter. 


2.3  The  Range  Gate  and  Transversal  Filter 

Analysis  of  the  transversal  filter  is  most  readily  carried  out  in  terms  of  the 
complex  envelopes  of  the  modulated  signals.  As  developed  in  section  3.3,  let 
r5  (t)  denote  the  complex  envelope  associated  with  r5(t).  The  range  gate  serves 
to  select  successive  pulse  returns  from  a  particular  range  cell.  Let  Ts  denote  the 
length  of  the  pulse  repitition  interval  (i.e.;  the  time  duration  between  consecutive 
pulses).  If  r5(t)  represents  the  current  pulse  return  from  a  particular  range  cell, 
then  the  (k-1)th  previous  return  from  the  same  range  cell  is  given  by  rs[t-(k-1  )TS]. 
The  corresponding  complex  envelope  is  denoted  by  r5  [t  —  (k-  1)Tg] . 

Assume  N  pulses  are  received  from  a  particular  range  cell.  The  transversal 
filter  sums  the  weighted  returns  to  produce  the  complex  envelope  output 

N 

r6(t)  =  X  wk*  -  r5 [t  —  (k-  1)T§]  .  (2.9) 

k  =  i 


wk  represents  the  kth  weight  where  the  tilde  denotes  a  complex  quantity  and  the 
asterisk  denotes  the  complex  conjugate.  A  block  diagram  of  the  transversal  filter 
is  presented  Fig. 2.6.  When  t  =  NTS,  the  output  is  given  by 

N 

’r6(NTs)  =  £  wk*  •  r5[(N-k+1)Ts]  .  (2.10) 

k  =  i 


Let 


"r5k  =  'r5[(N-k+1)TsJ  ;  k  =  1 ,2 . N  (2.11) 


and  define  the  vectors 


(w  ) 


T 


(2.12) 
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Figure  2.6  -  Block  Diagram  of  Transversal  Filter 


where  T  denotes  transpose.  Eq.(2.10)  can  then  be  expressed  as  the  inner  product 
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Chapter  3 

NONLINEAR  RECEIVER  ANALYSIS 


The  transmitted  radar  waveform  is  a  periodic  signal  composed  of  a  series  of 
sinusoidal  bursts,  whose  envelope  is  shown  in  Fig.  3.1.  Note  that 


Envelope  of  emitted  signal 


Figure  3. 1  -  Transmitted  Radar  Signal. 


E  =  pulse  amplitude 

fc  s  sinusoidal  carrier  frequency 

xe  =  pulse  width 

Ts  =  pulse  repitition  interval. 

An  analytic  expression  for  the  transmitted  signal  is 


E.cos(2jrfct) 


nTs+xe<t<  (n+1)Ts 


where  n 
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Radar 


3.1  Signal  at  Receiver  Input  (node  11 


components  of  v,  respectively.  The  carrier  frequency  of  the  received  signal  is 
(fc+fD)  where  the  Doppler  frequency  shift  due  to  the  radial  motion  of  the  target  is 
given  by  [4] 


2v  f 
r  c 


In  Eq.  (3.2)  c  denotes  the  speed  of  light  which  is  the  velocity  of  propagation  of  the 
electromagnetic  wave.  The  doppler  shift  plays  an  important  role  in  the  detection  of 
weak  signals  embedded  in  strong  clutter  returns. 

Assuming  that  the  delay  between  a  transmitted  and  received  pulse  is  known, 
there  is  no  loss  in  generality  by  omitting  the  delay  and  expressing  the  received 
signal  at  the  input  to  the  RF  ampifier  as 


^  (t)  =  as(t) cos [2jt(fc  +  fD)t +  6sJ +ad(t) cost 27ifct  +  0dJ  .  (3.3) 


Consider  a  received  signal  which  is  reflected  from  a  moving  target  with  ve¬ 
locity  vector  v,  as  illustrated  in  Fig.  3.2.  Let  vr  and  vn  denote  the  radial  and  normal 


Figure  3.2  -  Signal  Received  from  Moving  Target. 


li 


The  envelope  and  phase  of  the  desired  signal  are  denoted  by  a^t)  and  9S  while 
those  of  the  additive  disturbance  are  denoted  by  a<j(t)  and  0d .  In  terms  of  complex 
signals,  r-|  (t)  can  be  written  as 


(t)  =  Re{as(t)e 


j[2H(fc  + 


D,t  +  es] 


}  +  Re  {ad(t)e 


i[2ltfct  +  9d] 


Re  I  {as(t)e 


i<2It,D,  +  es)  jed, 

+  ad(t)e  }  e 


=  Re  {r1  (t)  e 


j  (27tf_t) 


} 


(3.4) 


where  Re{.}  denotes  the  real  part  operator.  The  complex  envelopes  of  the  desired 
signal,  disturbance,  and  received  signal,  respectively,  are 


s(t)  =  as(t)e 


d(t)  =  ad(t)e 


j<27tfDt  +  es) 

J0. 


(3.5) 


r-i  (t)  =  s(t)  +d(t). 


The  disturbance  consists  of  clutter  plus  thermal  noise.  Hence, 


d  (t)  =  c(t)  +  n  (t) 


(3.6) 


where  c  (t)  and  n  (t) ,  respectively,  denote  the  complex  envelopes  of  the  clutter 
and  noise,  s  (t) ,  c  (t) ,  and  n  (t)  are  each  modeled  as  statistically  independent 
zero  mean  stationary  complex  Gaussian  random  processes.  Because  a  sum  of 
jointly  Gaussian  random  processes  is  Gaussian,  d  (t)  and  r1  (t)  are  also  zero 
mean  stationary  complex  Gaussian  processes.  Hence,  at  a  specified  instant  of 
time,  their  magnitudes  and  phases  are  Rayleigh  and  uniformly  distributed  random 
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variables,  respectively.  In  addition,  each  phase  random  variable  is  statistically  in¬ 
dependent  of  the  corresponding  magnitude  random  variable.  r,(t)  is  assumed  to 
be  a  bandpass  signal  with  a  spectrum  as  shown  in  Fig.  3.3. 


Spectrum  of  rt{t) 


Figure  3.3  -  Power  Spectral  Density  of  rUt). 

3.2  Signal  at  RF  Amplifier  Output  (node  3) 

As  shown  in  Fig.  2.2,  the  RF  amplifier  is  modeled  as  a  zero-memory  nonlin¬ 
ear  amplifier  in  cascade  with  an  ideal  bandpass  filter.  The  output  of  the  nonlinear 
amplifier  is  determined  first. 


3.2.1  Output  of  Nonlinear  Amplifier  (node  2) 

As  given  in  Eq.(2.1),  the  input-output  characteristic  of  the  nonlinear  amplifier 
is 


r2  (t) 


(3.7) 


Since  the  nonlinearity  is  continuous  and  single-valued  and  has  all  derivatives  on 
an  interval  including  ^  =  0,  it  can  be  expanded  in  a  Mac  Laurin’s  series  at  r1  =  0. 
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Let  the  nonlinearity  be  denoted  by 


r2  =  9  0V 


The  Mac  Laurin’s  series  is  then  given  by 


r2  =  2  vrPi 


P  =  o 


where  the  p,h  coefficient  is 


_  1  dP 


zl  1  -3-5...(p-4)  (p - 2) 

o2  P! 

g 


\  p=i 


;  p=3,5,7,..  (3.10) 


;  p=0,2,4,6,.. 


Letting  p  =  2k+1 ,  it  follows  that 


~  f-1  f  (2k- 1)!!  2k  +  1 
1  A  o2  (2k+1)!  ' 

K  =  1  n 


(3.11) 


where  the  double  factorial  is  defined  to  be 


(2k -  1 ) !!  =  1  -3-5 . (2k-3)  (2k-  1 ) 


(3.12) 
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Eq.(3.11)  is  now  expressed  in  terms  of  the  complex  envelope  of  r-j  (t).  Recall 


j  A 

ri(t)  =  Re  {r,  (t)e  c  } 


j27lf  t  _  .  -j27Cf  t 

Mt)e  c+r1  (t)  e  c 
2 


(3.13) 


where  the  asterisk  denotes  complex  conjugate.  It  follows  that 


,2k  + 1 


1  [.  j2*U  .  *  -j27tf  t 

- r.e  +  r<  e 

,2k +  i  I  1  1 


2k+ 1 


(3.14) 


Application  of  the  binomial  expansion  to  Eq.(3.14)  results  in 


,2k +  i 


2k  +  i 

— C—  Y  Cm+1 
02k  +  i  ^  m 

4  m  =  o 


j27tfji2k  +  1  m  r  *  -j2^f  tim 
r,e  0  f,  e 


where  the  binomial  coefficient  is 


(3.15) 


,2k +  i  =  (2k+1)! 

'm  m!  •  (2k  +  1  -  m) ! 


(3.16) 


Combining  the  exponentials,  Eq.{3.15)  becomes 


2k  + 1 


r?k+1  =  -j— •  y  C 

1  n2k+1  " 


m  =  0 


2k+i  r-  ,2k+i-m 
m  '  lriJ 

*  m  j27ifrt(2k+i  -2m> 

•  (r,  ]  e  c  .  (3.17) 
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It  is  seen  that  raising  r-|  to  the  power  (2k+1 )  creates  spectral  zones  located  at  the 
frequencies  (2k+1-2m)fc  where  m  is  an  integer  ranging  from  0  to  (2k+1).  The 
spectral  zones  are  referred  to  as  even  or  odd,  respectively,  depending  upon 
whether  the  factor  multiplying  fc  is  even  or  odd. 


The  following  observations  are  now  made  with  respect  to  the  various  spec¬ 
tral  zones  : 


1 )  As  m  varies  from  0  to  (2k+1 ),  the  locations  of  the  spectral  zones 
range  from  (2k+1  )fc  to  -(2k+1  )fc. 

2)  All  of  the  spectral  zones  are  odd  because  (2k+1-2m)  is  always  an 


odd  integer. 

3)  An  odd  spectral  zone  occurs  at  the  positive  frequency  (2p+1  )fc 

when 


when 


(2k+  1  -2m)  =  (2p+1) 

m  =  k-p  .  (3.18) 

4)  An  odd  spectral  zone  occurs  at  the  negative  frequency  -(2p+l)fc 


(2k  +  1  -2m)  =  -(2p+  1) 

m  =  k  +  p  +  1  .  (3.19) 


5)  Zones  located  at  positive  frequencies  result  when  (2k+1-2m)  >  0. 
This  will  be  the  case  provided 


0<m£k  . 


(3.20) 


Substituting  Eq.(3.18)  into  Eq.(3.20)  and  simplfying  produces  the  inequality 

0  <  p  <  k  .  (3.21) 

6)  Zones  located  at  negative  frequencies  result  when  (2k+1-2m)  <0. 
This  will  be  the  case  provided 

k  +  1  <,  m  <2k+  1  .  (3.22) 
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Substituting  Eq.(3.19)  into  Eq.(3.2 2)  and  simplifying  also  produces  the  inequality 

O^psk  .  (3.23) 

Using  the  above  observations,  Eq.(3.17)  is  readily  expressed  as  two  sum¬ 
mations  where  the  first  sum  is  over  negative  frequencies  and  the  second  sum  is 
over  positive  frequencies.  In  particular, 


r2k+ 1 

ri 


r2k  +  i 
uk+p+i 


k  +  p+  i 


-j27I(2p  +  1)f  t 

e  c 


k 

*10 


p  -  ° 


2k+ 1 
k-p 


fr1]k  +  P  +  1 


] 


k-p 


j27C(2p+1)f 


e 


'} 


ib(  l 

*  n  -  i 


2k  +  i|  "  k  +  p  +  i 

p  =  0 


B-j*,(*p+„»ct+  ^  ,  .[?i]k-p 


P  =  o 


[f,]2p+1  ■  e 


j2]t(2p 


+  1),c'  J 


-Pi 

*  n  —  n 


«2k  + 1 

2  *p  =  0 


:  ,2(k-p)  [- -]2p  +  1 


-j2*<2P+1)fct  k  k  +  ,  2(k_p) 

6  +  L  uk-p  'lril 

p  =  o 


-  2D  +  1  j2K(2p+1)fct 
[r1]^p  1  •  e  c 


} 


(3.24) 
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Evaluation  of  the  binomial  coefficients  results  in 


*2k  +  i 


(2k+  1)! 


(2k+  1)1 


k  +  p  +  i  (k  +  p+1)!  •  (2k +  1  —  k  —  p  —  1 ) !  (k  +  p  +  1)!  ■  (k-p)! 


»2k+  i 


(2k  +  1 ) ! 


'k-p  (k-p)!  •  (2k  +  1  -k+p)! 


(2k  + 1 ) ! 

(k-p)!  •  (k  +  p  +  1)! 


(3.25) 


It  is  concluded  that 


s2k  + 1  _  r2k+ 1 

'k  +  p  +  i  ~  °k-p 


(3.26) 


Eq.(3.24),  therefore,  simplifies  to 


,2k +  i 


^r{i^k-VN2<k'P)[rri*]2pt1 

vp  =  0 


e 


_j2it  (2p  + 1 )  ft 


c  +  [f1]2p+’ 


j27E(2p 


+1)Vj  J 


2k+i  i-i2(k-p> 


1  .  y  c2k+1  •  If 

2k  ^  Gk-P  I  1 
2  p  =  0 


i27tV 


2p  +  1 


(3.27) 


Substituting  Eq.(3.27)  into  Eq.(3.11),  the  expression  for  the  nonlinear  amplifier 
output  becomes 


oo  /  .k  u  r>2k+ 1 

.  .  V  f=l)  V  IHzP.  I?  l2(k-P) 

1  ^  I  a2  (2k  + 1)!  ^  „2k  I  11 


p  =  o 


Jr  ]2nf  tfP  +  1 

Re  r^  c 


(3.28) 
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Separating  out  the  spectral  zones  located  at  ±fc,  the  above  expression  can  be  re¬ 
written  as 


r„  = 


2L 


I 


2  no2 
9 


oo 

f"1) 

_2 

1 

k=  1 

l  gJ 

K  n2k  + 1 

°k  (2k -  1 ) !!  |2k 


,2k  (2k  +  1)! 


i 


Re  {r.,e  c} 


2L 


I 


2jio‘ 


g  L 


00  ^_i  ^ 


s 

k  =  i 


(2k-  1 ) !! 


k  C?k+1 

S  k-p 


.2  (k  —  p> 


(2k+1)!  ^  02k 

p=  i  * 


Re 


r.  )2itfcr 

2p  +  v 

1  1 

r,e 

J 

[J 

(3.29) 


Let  Ar2  denote  the  terms  in  r2  corresponding  to  the  spectral  zones  other  than 
those  centered  at  ±fc.  In  addition,  define  the  complex  envelope 


r2 


2L 

feno2 

V  g 


OO 

I'-l'l 

1+ 1 

II 

_ 1 

t°gJ 

(2k-  1 ) !!  2k  - 

(2k+  1 ) !  Iril  f 


(3.30) 


It  is  then  possible  to  express  the  nonlinear  amplifier  output  as 


' 

?2 

•  e  l 

f 

. 

c 

(3.31) 


The  spectrum  of  r2(t)  is  illustrated  in  Fig. 3.4.  The  spectral  zones  at  ±3fc,  ±5fc>... 
arise  from  terms  in  Ar2.  The  spectral  zones  at  ±fc  are  due  to  the  first  term  in 
Eq.(3.31). 
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Spectrum  of  r2(t) 


Figure  3.4  -  Power  Spectral  Density  of  r0  ft). 

3.2.2  Output  of  Bandpass  Filter  (node  3) 

The  bandpass  filter  is  assumed  to  have  the  ideal  transfer  function  whose 
magnitude  is  shown  in  Fig. 2.4.  The  bandwidth  Brf  is  assumed  to  be  narrow 
enough  such  that  only  the  spectral  zones  at  ±fc  are  passed  with  negligible  distor¬ 
tion.  The  remaining  spectral  zones  are  assumed  to  be  perfectly  attenuated.  It  fol¬ 
lows  that  the  output  of  the  RF  amplifier  is  given  by 


r3(t) 


=  ReJ 


* 

i^c'l 

.  ?2 

•  ®  > 

'c  J 

(3.32) 


where  the  complex  envelope 
shown  in  Fig. 3.5. 


is  defined  in  Eq.(3.30).  The  spectrum  of  r3(t)  is 


3.3  Signal  at  Mixer  Output  (node  5) 

As  discussed  in  Section  2.2,  the  mixer  consists  of  the  multiplier,  local  oscilla¬ 
tor,  and  lowpass  filter.  The  input  to  the  lowpass  filter  is 

r4  (t)  =  r3(t)  •  C°s(27cfct)  •  (3.33) 
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In  addition, 

n  r  j2;if  t  -j2jtf  r 

cos(27if  t)  =  -  •  e  c  +e  c  •  (3.35) 

c  2  L 


It  follows  that 


r4 (t)  =  X  r2f  *e  +r2  f  e 
L  'c  'c 


i  i  - .  -  *  i  r  -i27cfcr 

1 '  1  n  c  +r9  -e  c  •  e  c  +e  c 


1  r.  .  *  1  1  f.  J27t  (2f  )t  -J27C  (2f  )t 

i  •  r2  +  r2  +  —  ■  r2  •  e  +  r2  •  e 

4  L  L  4  L  2  L 

L  c  cj  L  c  c 


1  .  if-  j2K(2f  )t 

2  •  Re {r2  }  g  ’  Re V 2 f  *® 

Tc  rc 


(3.36) 
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Observe  that  the  spectral  zones  of  r4(t)  are  located  at  dc  and  ±2fc.  The  spectrum 
of  r4(t)  is  shown  in  Fig. 3.6. 


Spectrum  of  r4(i) 


The  lowpass  filter  is  assumed  to  have  the  ideal  transfer  function  whose  mag¬ 
nitude  is  shown  in  Fig.2.5.  The  bandwidth  BL  is  assumed  to  be  wide  enough  to 
pass  with  negligible  distortion  the  spectral  zone  at  dc  but  small  enough  such  that 
the  spectral  zones  at  ±2fc  are  perfectly  attenuated.  It  follows  that  the  signal  at  the 
mixer  output  is 

r5  (t)  =-■  Re  {r2  }  .  (3.37) 


Define  the  complex  signal 


(3.38) 


where  r2 


is  defined  in  Eq.(3.30).  Then 


r5(t)  =  Re (r5(t) }  - 


(3.39) 
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The  spectrum  of  r5(t)  is  illustrated  in  Fig.3.7.  Note  that  r5(t)  is  a  base-band  signal 
which  serves  as  the  input  to  the  range  gate  and  transversal  filter. 


Spectrum  of  r5(t) 


Figure  3. 7  -  Power  Spectral  Density  of  fc(t). 
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Chapter  4 


RELATIONS  BETWEEN  AUTOCORRELATION 
FUNCTIONS  AT  VARIOUS  POINTS  OF  THE 

RECEIVER 


In  this  chapter  the  general  form  of  the  complex  autocorrelation  function  at 
the  input  to  the  transversal  filter  (node  5  of  Fig. 2.1 )  is  derived  as  a  function  of  the 
complex  autocorrelation  function  at  the  output  of  the  zero-memory  nonlinear  am¬ 
plifier  (node  2  of  Fig. 2.2). 

4.1  Autocorrelation  Function  at  the  Input  to  the  RF  Amplifier  (node  1  of 
Fie.2.1) 

The  autocorrelation  function  of  the  signal  r-j(t)  at  the  input  to  the  RF  amplifi¬ 
er  is  defined  as 


R^tt)  =  E  [r,  (t)  •  r,  (t  —  x)  ]  . 


(4.1) 


In  terms  of  the  complex  envelope  r1  (t)  introduced  in  Eq.(3.4),  the  autocorre¬ 
lation  function  becomes 


F^at)  =  E 


j27Cf  t  _  j27tf  (t-T) 

Re {r.,  (t)e  c  }  •  Re  (r,  (t-x)e  c  } 


=  E 


j2ref  t  .  *  -j2Jtf  t 
^  (t) e  c+r1  (t) e  c 


j27tf  (t-T)  -j27tf  (t-T) 

r1  (t-x)e  c  +  r1  (t  —  x)  e 
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Rl  (t,X)  = 


2 

2 


j27lf  (2t-X)  _  *  -j27lf  (2t  -  T) 

n  (t)  r,  (t-x)e  c  +r1  (t)^  (t-x)e 

2 


j2Jtf  T  _  .  ,  -j27tf  T  "I 

(t-x)  ■  e  +r1  (t) -r,  (t-x)  e 
2 


Interchange  of  the  real  part  and  expectation  operators  results  in 


1  .  -  j2JCf  (2t-X) 

FMt.x)  =  5-Re{E[r1(t)  rt  (t - x) ]  e  c  } 


\2rtcT. 


+  2  •  Re(E[ri  (t)  •  ri  (t - x) ]  •  e  }  -  (4.2) 


Assuming  stationarity  of  r-|(t)  within  the  gating  window,  it  is  shown  in  Appendix  A 
that 


E[ri(t)  -  ?!  (t  -  x)  ]  =  0 


(4.3) 


In  addition,  stationarity  implies  that  E  [  r1  (t)  •  r1  (t  -  x)  ]  is  a  function  of  the  differ¬ 
ence  between  the  arguments.  Hence,  we  define 


Ri  (x)  =  E [?!  (t)  •?/  (t-x)]  • 


It  follows  that 


R1  (t,  x)  =  R,(x)  =  ^  ' 


Re 


Ri  (x)  ■  e 


j27tfcx- 


(4.4) 


(4.5) 
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4.2  Autocorrelation  Function  at  the  Output  of  the  Zero-memor 


Nonlinear  Amplifier  (node  2  of  Fis.2.2 


4.2.1  Transformation  of  the  probability  Density  Function  Due  to  the 


Nonlinearity  f51 


As  pointed  out  in  section  2.1.1,  the  nonlinear  amplifier  is  modeled  as  a  soft 
nonlinearity  with  saturation  levels  at  ±L.  In  particular,  as  given  by  Eq.(2.1),  the  in¬ 
put-output  characteristic  is  described  by 


r’(t)  /  2  \ 

r2(t)  =  |  exp  dZ 

°  t 2<v 


Let  the  probability  density  function  of  r^t)  and  r2(t)  be  denoted  by  p-|(r-,)  and 
p2(r2),  respectively.  These  two  densities  are  related  by  the  transformation 


From  Eq.(2.2), 


P2(r2)  =  Pi  (ri)  ' 


Because  the  slope  is  nonnegative, 
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Iri  =  g  (r2) 


12710 ‘  2G2 

_ g ,e  g 

2L  6 


(4.7) 


P2(r2)  = 


rKOn  2°2 

~2\~  '  Pi  (ri )  '  e  9  -1 

2L  r,  =  g  (r2> 


(4.9) 
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Consider  the  case  for  which  r,(t)  is  a  zero-mean  Gaussian  random  variable 
with  variance  c^2.  Then 


Pi  (ri) 


It  follows  that 


(4.10) 


P2(r2)  = 


A  r?  1 

J2nGn 

2a2  2G2 

*  g  p 

2L  ® 

k  g  \) 

ri=g  (r2) 


ag/0i 

2L 


=  g_1  (r2) 


(4.11) 


Recall  from  Eq.(2.7)  that  the  nonlinearity  can  be  approximated  by  a  linear  charac¬ 
teristic  over  the  range  (-ag,+ag),  Also,  note  that  o-i  is  the  rms  value  of  the  input 
r^t).  Hence,  the  ratio 


(4.12) 


is  a  measure  of  how  weakly  the  nonlinearity  is  being  driven  by  the  input  signal. 
The  larger  is  the  value  of  a,  the  more  likely  it  is  that  the  input  will  be  in  the  linear 
region  of  the  nonlinearity.  In  terms  of  a,  the  density  function  at  the  output  of  the 
nonlinearity  is  given  by 


a 


a 


2a. 


W*)]2 


P2<r2)  =  2L  0 


g 


(4.13) 
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The  nature  of  this  density  function  is  now  discussed  for  various  values  of  a. 
When  a  =  1, 


P2(r2) 


;  -L  <  r2  <  L 
;  otherwise  . 


(4.14) 


A  plot  of  Eq.(4.1 4)  is  shown  in  Fig.4.1 .  We  conclude  that  the  nonlinear  amplifier 


pm 


Figure  4. 1  -  Plot  of  Pofro)  when  a  =  1. 

output  is  uniformly  distributed  when  cr1  =  ag  .  For  a  >  1 ,  the  exponent  in  p2(r2)  is 
negative. From  Fig. 2. 3,  it  is  obvious  that 

g  1  (-L)  =  anc*  g  1  (L)  =  oo 

Consequently,  P2(r2)  is  a  bell-shaped  curve  whose  value  asymptotes  to  zero  for 
r2  =  ±L.  For  a  »  1 ,  the  nonlinearity  is  exercised  primarily  in  its  linear  region,  it  is 
not  suprising,  therefore,  that  the  corresponding  density  function  is  closely  Gauss¬ 
ian.  However,  p2(r2)  cannot  be  strictly  Gaussian  because  r2  is  constrained  to  the 
interval  (-L,L)  whereas  a  Gaussian  random  variable  can  assume  values  over  the 
infinite  interval  (-«,+<«).  For  a  <  1 ,  the  exponent  in  p2(r2)  is  positive.  It  follows  that 
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the  density  function  becomes  infinite  for  r2  *  ±L.  In  the  limit,  as  a  — » 0,  the  density 
function  approaches  two  impulses  located  at  r2  =  ±L,  as  shown  in  Fig. 4.2.  This  is 


P2(r2> 


Figure  4.2  -  Plot  of  Do(ro)  in  the  limit  as  a  ->  0. 

consistent  with  the  fact  that  the  input-output  characteristic  approaches  that  of  the 
hard  limiter,  illustrated  in  Fig.4.3,  as  ag  — »  0.  Even  if  ag  *  0,  the  nonlinearity  is 
driven  into  saturation  most  of  the  time  when  a-j  »  ag  or,  equivalently,  a  «  1  The 
discrete  density  function  of  Fig. 4.2  is  then  applicable. Sketches  of  p2(r2)  for  vari¬ 
ous  values  of  a  are  shown  in  Fig.4.4.  From  Fig. 2.3,  note  that 

g-'  (0)  =  0  ■ 


Consequently,  the  curves  in  Fig.4.4  obey  the  property  that 

P2<°>  ^  It  ’ 


(4.15) 


Also,  as  pointed  out  previously, 


if  a>  1 
if  a  <  1  . 


(4.16) 
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Fiaure  4.3  -  Input-output  characteristic  of  hard  limiter. 


Fiaure  4.4  -  Sketches  of  Pofro)  for  various  values  of  a. 


4.2.2  Transformation  of  the  Correlation  Function  Due  to  the 


Nonlinearity  f51 


Define  the  Fourier  transform  of  the  input-output  characteristic  to  be 


G(ju)  =  J  g [r1  (t)  1  - e“^Uri <t>dr1  (t)  . 


(4.17) 


Assuming  stationarity,  the  output  correlation  function  is  given  by 


R2(t)  =  E[r2(t)  r2  (t - x) ] 


(4.18) 


The  output  voltages  can  be  expressed  in  terms  of  the  Fourier  transform  of  the  in¬ 
put-output  characteristic  according  to  the  relations 


r2(t)  =  g[rt  (t)  ]  =  J-  j  G(ju)  •  e‘ 


jur,  (t) 


(4.19) 


oo 

1  r  ^  .  jvr,  (t-t) 
r2  (t  —  x)  =  g  [r1  (t  —  x)  ]  =  J  G  (jv)  eJ  < 


(4.20) 


It  follows  that 


R2  (X)  „  e[ -L  |  J  G  (ju,  .  G  (jv)  ■  eiur’ (,)  •  j"' " '  T)  dudv 

_  —OO—OO 


JL  J  J  G(ju) -G(jv)  E[ejuri(t)  ejvri(t  T)]dudv  -(4.21) 


The  quantity, 


.  cr  jur, (t)  jvr, (t-x)i 
M(ju,jv)  =  E[e  e  j 


(4.22) 


is  recognized  as  the  joint  characteristic  function  of  the  random  variables  r,  (t)  and 
r-i(t-T).  Hence, 


R2  (*)  =  -^5  f  f  G  (ju)  •  G  (jv)  •  M  (ju,  jv)  dudv 
4  nz  L  L 


(4.23) 


To  evaluate  GQ'u),  use  is  made  of  the  differentiation  property  of  the  Fourier  trans¬ 
form.  To  begin  with,  we  have  the  Fourier  transform  pair 

n  (t) 

exp 
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Applying  the  differentiation  property  to  Eq.(4.27),  there  results 


(t)) 
dr1  (t) 


2L  f  ri  W  1 

—  exp - —  ju  •  G  (ju) 


(4.28) 


However,  it  is  well  known  that  the  Fourier  transform  of  a  Gaussian  waveform  is  a 
Gaussian  waveform,  as  indicated  below: 


2  2 

cgu 


fr1(t)l  - 2  2 

exp  J2Jt<  e 


(4.29) 


It  follows  that 


o2  u2 
9 


ju  •  G  (ju)  = 


2 no2  e 


(4.30) 


Thus, 


G(iu)  =  |e  2 


(4.31) 


Substituting  Eq.(4.31)  into  the  expression  for  R2(t),  we  have 


(2d27  7-i  -fiuw, 


’l  r.l  ..2 


[U  +v  ]  -R,  (t)uv 


dudv 


(4.32) 
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Combining  terms  in  the  exponent  of  the  integrand, 


(2L) 


R2(x)  =  -5-4-  f  f  — 

2  4tc2  J  J  uv 


-1  l 


[u2  +  v2]  -  R,  (t)uv 


e 


dudv 


(4.33) 


To  obtain  an  integrand  without  the  factor  ,  we  differentiate  R2(t)  with  respect 
to  Ri(t).  Specifically, 


dR2(t)  (2L)‘ 


ra2+a4 


dR1  (t)  4^2 


J 1 


—  00—00 


•  [u2  +  v2]  -R^tjuv 


e 


dudv  -  (4.34) 


At  this  point,  it  is  convenient  to  make  the  change  of  variables 


x  = 


2  2 
G  +07 

9  1 


V-T- 


u  and 


V= 


2  2 
Ot  +  G* 


(4.35) 


This  results  in 


dR2(x)  _  (2L)2  2 

_____ 


2  2  2R,  (X) 

c  oo  -<X  +y)-Z^-Z2 

G  +  O 


xy 


J  S 


4n2  o2  +  o2 

g  1  -oo-oo 


e 


9  dxdy  •  (4-36) 


To  further  simplify  the  integrand,  introduce  the  variable  v  and  y  such  that 


x  =  h(u,y)  =  x>- 


y  =  k(u,  y)  = 


i 


(4.37) 


J 1  -  a2 
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where 


a  = 


R1  (t) 


2  2 
C  +CT* 


g 


Observe  that 


o  =  x  +  ay 


y  =  J 1  -a2y  - 


The  Jacobian  of  this  transformation  is  given  by 


J  (v,  Y) 


dh  dh 

[i 

du  dy 

a/i  -  a2 

dk  dk 

du  dy 

[h-a2\ 

(4.38) 


(4.39) 


(4.40) 


Hence, 


dxdy  = 


Ji-? 


dodY  . 


(4.41) 


Utilizing  the  transformation  of  Eq.(4.37),  the  exponent  of  the  integrand  in 
Eq.(4.36)  becomes 


-[x2  +  y2  +  2axy]  =  - 


2a  -v.  a 


J^a2 


1  -  a 


2  2a 

;t  +  -r—vy 

VI  -  a2 


1 


1  -  a‘  1  -  a‘ 
-[^h-y2]  - 


2  a2  + 1  -  2a2 

U  + - 2 

1  -  a2 


Y2 


(4.42) 


Thus,  the  expression  for 


dR2(x) 

dR^t) 


becomes 
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dR2<T>  (2L)2  2  1  7  T.-l^  +  T2].  . 

hr-t;v  =  - o - 2 - 2  '  -7=5  I  I  e  dudY  •  (4-43) 

dRi^  4k2  o2+o2  /i  _a2  J  J 

g  I  \  1  Q  —oo—oo 


The  double  integral  is  readily  evaluated  to  yield 


J  J  e  ^  +^dudy  =  Jt  . 


(4.44) 


Consequently, 


dR2(t) 

dRTCO 


2  2 
Cg  +  °1 


R-l  (X) 

a  +of 
a  i 


(4.45) 


Recalling  that 


d[sin  nu] 


1  du 

L  2  dx 
/I  -  u 


(4.46) 


it  follows  that 


RP(x)  = 


(2L)2  .  ^(t) 

-sin  1  -J — 5  +C 
271  o2  +  o2 

L  g  iJ 


(4.47) 


where  C  is  a  constant  of  integration.  To  evaluate  C,  consider  the  special  case  for 
which 


a2  =  o2  and  x  =  0 
9  1 


(4.48) 


Then 
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R,(0)  = 


R,(0)  = 


(2D 


(4.49) 


Because  ag2  =  oi2,  a  =  1 .  It  follows  that  p2(r2) is  the  uniform  density  function  giv¬ 
en  by  Eq.(4.14).  As  a  result, 


B2(0)  =E[ri]  = 


(4.50) 


By  comparison  of  Eqs.(4.49)  and  (4.50),  it  is  concluded  that  C  =  0.  The  output  cor¬ 
relation  function,  therefore,  is  given  by 


r2w  =  ^ 


(2L)2  r  Ri (x) 


2  ,  2 
°g  '  J 


(4.51) 


Recalling  that 


R1  anH 

P,  (x)  =  — 2 —  and 


o 

_  g 


R2(x)  can  also  be  expressed  as 


„  (2L)2  -TPi  <*>" 

R2(X>  =  pZ  "  Sin  - 2  ■ 

2 *  ^+a2 


(4.52) 
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4.3  Autocorrelation  Function  at  the  Output  of  the  RF  Amplifier  (node  3 
nf  Fig.2.1) 


From  Eq.(3.32)  the  RF  amplifier  output  is 


r,(t)  =  Re 


h  (t) 


e 


j2itfct'| 


(4.53) 


where  r2  (t)  is  given  by  Eq.(3.30).  By  definition,  the  autocorrelation  function  of 
fc 

r3(t)  is 


R3(t,T)  =  E[r3(t)  r3  (t  -  x)  ] 


(4.54) 


Following  a  development  analogous  to  that  which  resulted  in  Eq.(4.2), 


R3(t,x)  =  2Re 


r2(t) 


r2(t-x) 


f 


j27Cf_  (2t-X) 


e 


1 

2 


- 

_  * 

- 

j27tf  XI 

(E 

r2  (t) 

r2  (t-x) 
c 

fc. 

•  e 

(4.55) 


In  Appendix  B  it  is  shown  that 


Mt) 


r2  (t  -  x) 


f. 


=  0 


(4.56) 


and  E 


r2(t)  f  •  r2*  (t  —  x) 
c 


'cJ 


is  a  function  only  of  x.  Define 


R2(t) 


=  E 


r2(t) 


-  * 


r2  (t-x) 


'  '2 


Cj 


(4.57) 


38 


The  expression  for  R3(t,t)  then  simplifies  to 


1  -  i2*fcx 

R3(t,t)  =  R3(x)  =  «  •  Re{R2(T)  e  c  } 

*  Tc 


(4.58) 


4.4  Autocorrelation  Function  at  the  Mixer  Output  (node  5  of  Fig.2.1 


As  pointed  out  in  section  3.3,  the  signal  at  the  mixer  output  is  a  baseband 
signal  given  by 


r5(t)  =  Re  {r5  (t) } 


(4.59) 


where 


h  (*)  =  ,•  h  «> 


(4.60) 


The  autocorrelation  function  of  r5(t)  is 


r  r2(,,|  .  r2  (t  -  t)  | 

R5(t,T)  =  E[r5(t)  r5(t-t)]  =E  Re  — ^  •  Re  - ^ - - 


lfc  ‘  lfc. 


=  ?2  (t)  +?2  (t)  ?2  (t  -  X)  +?2  (t-T)  } 


fr2(t)  r2 (t -x)  +  r2  (t)  r2  (t-x) 

1CJ  *c  fc  Tc  *c 

8  ]  2 


r2(t)  r2  (t-x)  +r2  (t)  r2(t-x) 

*C  *C  *C  *C 
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=  qE {Re  r2 (t)  r2 (t  —  x)  +Re  r2 (t)  r2  (t-t)  } 


lfc  I'c 


f'c  lfc. 


Re {E  ?2(t)  r2 (t - T)  +  E  r2(t)  r2  (t-t)  } 

c  cj  L  c  c_ 


(4.61) 


As  noted  in  section  4.3, 


E  r2(t)  •  r2(t-x)  =  0 

Tc  'c_ 


(4.62) 


and  E  r2  (t)  •  r2  (t  - 1)  .  is  a  function  only  of  x.  With  reference  to  Eq.(4.57), 

fc  fcJ 

define 

R2  (^)  |f 

R5(t)  =  Ie^jO)  ?2  (t-t)  1  =  - J- -  •  <4  63) 

L  Tc  'cJ 


It  follows  that 


R5(t,x)  =  R5(t)  = 


gRe{R5(t)}  =  gRe{R2(t)|f  } 


(4.64) 


Equation  (4.63)  gives  a  direct  relationship  between  the  complex  autocorrelation 
functions  at  the  outputs  of  the  mixer  and  zero-memory  nonlinear  amplifier. 
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Chapter  5 

AUTOCORRELATION  FUNCTION  AT  THE 
TRANSVERSAL  FILTER  INPUT  FOR  THE 
GAUSSIAN  CASE 

The  results  of  Chapter  4  are  used  in  this  chapter  to  derive  an  expression  for 
the  complex  autocorrelation  function  at  the  input  to  the  transversal  filter.  If  s(t),  c(t) 
and  n(t)  denote  the  desired  signal,  clutter  and  noise,  respectively,  the  received 
signal  is  given  by 


r,  (t)  =  s(t)  +c(t)  +n(t)  .  (5.1) 

The  sum  of  clutter  and  noise  is  referred  to  as  the  disturbance,  d(t).  Thus, 

d(t)  =  c(t)  +  n(t)  .  (5.2) 

The  desired  signal,  clutter  and  noise  are  assumed  to  be  statistically  independent 
zero-mean  Gaussian  random  processes.  As  a  result,  n(t)  is  also  a  zero-mean 
Gaussian  random  process. 

5.1  Autocorrelation  Function  at  Receiver  Input 

Assuming  stationarity,  the  autocorrelation  function  of  the  received  signal  is 

R1  (x)  =  E  [r,  (t)  •  r1  (t- x) ] 

=  E  {  (s  (t)  +c(t)  +n(t)]  •  [s  (t  —  x)  +  c(t-x)  +  n(t-x)]}  -(5.3) 

Because  of  statistical  independence  and  zero-means,  the  expectations  involving 
cross  terms  are  zero.  Therefore, 

R-j  (t)  =  Rs  <x)  +  Rc  (x)  +  Rn  (t)  (5.4) 
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where 


Rs(x)  =  E [s (t)  •  s(t-x)] 


Rc(t)  =  E  [c  (t)  •  c  (t  —  x)  ] 


Rn(t)  =  E  [n  (t)  •  n  (t  —  x)  ]  . 

The  desired  signal,  clutter  and  noise  can  be  expressed  in  terms  of  their  complex 
envelopes  as 

j27tf  t 

s (t)  =  Re{s(t)  e  c  } 


c  (t)  =  Re  {c (t)  •  e  c}  (5-6) 


n (t)  =  Re{n(t)  e  u  }  . 


Following  the  same  derivation  which  was  used  to  obtain  Eq.(4.5)  from  Eq.(4.1), 
the  various  autocorrelation  functions  are  related  to  their  complex  autocorrelation 
functions  by 

1  -  j27tf  x 

Rs(t)  =  ^  Re{Rs(x)  e  c  } 

1  .  j27Cf  X 

Rc(x)  =^Re{Rc(x)e  c}  (5.7) 

1  -  J27Cf  X 

Rn(x)  =  ^  •  Re {Rn (x)  •  e  c  } 
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where 


Rs(x)  =  E  [s (t)  •  s  (t - x) ] 

Rc(x)  =  E[c(t)  c  (t-x)  ]  (5.8) 

Rn(x)  =  E[n(t)  •  n*  (t-x)]  . 

Substitution  of  Eqs  (5.7)  into  (5.4)  results  in 

1  - 

R-i  (x)  =  I-RefR^x)  e  c  } 

1  j27lf  t 

=  ^■Re{[Rs(x)+Rc(x)+Rn(x)]  e  c}  -  (5-9) 

It  follows  that 

Ri(x)  =  Rs(x)+Rc(x)+Rn(x)  .  (5.10) 


From  Eq.(3.5)  the  complex  envelope  of  the  desired  signal  is 

j(27ifDt+es) 


S  (t)  =  as  (t)  •  e 


Thus,  the  complex  autocorrelation  function  of  s  (t)  is 


Rs(x)  =  E{a  (t)  e 


j(2nfDt  +  0  ) 


ag (t-x)  •  e 


-][2jcfD(t-x)+0s] 


=  E 


j27lfDX 

(as(t)  as(t-x) ] e  T 


(5.11) 


Assuming  the  fading  to  be  constant  during  the  interval  over  which  the  N 
pulses  from  a  given  range  cell  are  received,  as(t)  can  be  approximated  by  a  ran¬ 
dom  constant  A  which  is  Rayleigh  distributed.  Hence, 
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(5.12) 


E[as(t)  as(t-x)]  =  E  [A^] 


Use  of  Eq.(5.12)  in  Eq.(5.11)  results  in 


Re(x)  =  E  [A  ]  e 


2,  J2n'DX 


(5.13) 


Let  as2  denote  the  variance  of  s{t).  From  Eqs  (5.7)  and  (5.13)  note  that 


Rs(°)  =  a2  =  gRe  (Rs(0) }  = 


E[A2] 


(5.14) 


Therefore, 


E  [A2]  =  2o2  . 


(5.15) 


Consequently,  the  complex  autocorrelation  function  of  s(t)  becomes 


P  j27cfDx 
Rs(x)=2o|-e 


(5.16) 


Substitution  of  Eq.(5.16)  into  Eq.(5.7)  results  in 


Rs(t)  =  2a|  - cos[27c(fc  +  fD)x] 


(5.17) 


As  shown  in  Fig. 5.1,  the  clutter  power  spectral  density  is  assumed  to  be 
Gaussian  shaped  around  the  carrier  frequency.  Consequently,  the  autocorrelation 
function  of  the  clutter  is  modeled  as 


Rc(t)  =o2  e  w  •  cos  (2nfcx) 


(5.18) 
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m 


Figure  5. 1  -  Gaussian  shaped  clutter  spectrum. 
The  corresponding  power  spectral  density  is 


Sc(f)  = 


{ 


f 


-2f7CO  ..(f-f 


•  2ko^  ■  e 
w 


w 


i -  -2r7to  (f+f_)T 

02JW.e  L  w  cJ 

c  \  w 


;f  >0 
;f  <0 


(5.19) 


Note  that 


o2  =  R  (0) 
c  c  ' 


(5.20) 


is  the  variance  of  c(t).  Also,  cw2  is  the  variance  of  the  Gaussian  function  used  to 
model  the  clutter  autocorrelation  function  envelope.  As  such,  it  is  a  measure  of  the 
width  of  the  clutter  autocorrelation  function.  The  larger  is  aw2,  the  wider  is  the  clut¬ 
ter  autocorrelation  function  and  the  narrower  is  the  clutter  power  spectral  density. 
Rc(t)  can  be  expressed  as 


Rc(x) 


20 

o2  •  e  w 
c 


j27tf  X 

•e  c 


(5.21) 
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Comparison  of  Eq.(5.21)  with  Eq.(5.7)  results  in 


Rc(t)  =  2o*-e  2°w  -  (5-22) 

The  additive  noise  is  assumed  to  be  bandlimited  white  noise.  Hence,  its 
power  spectral  density  is  a  constant  over  a  band  of  frequencies  of  width  W  cen¬ 
tered  at  ±fc,  as  seen  in  Fig. 5.2.  The  corresponding  autocorrelation  function  is 


Bn(t)  =  j  Sn  (f)  e 


=  o2 
n 


sin  (7tWx) 
7tWt 


cos  (27tfcx) 


(5.23) 
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By  comparing  Eqs(5.7)  and  (5.24),  it  follows  that 


Rn(t) 


O  2  Sin(TtWT) 
n '  reWx 


(5.25) 


Combining  the  expressions  for  Rs  (x),  Rc  (x)  and  Rn  (x),  we  obtain 


R-i  (x)  =  2of  •  e 


j2JCfnX 


20 


+  2  o“  •  e 


w  .  „  2  sin  (icWx) 


2o 


7tWx 


(5.26) 


5.2  Autocorrelation  Function  at  the  RF  Amplifier  Output 

The  autocorrelation  function  at  the  nonlinear  amplifier  output  (node  2  of 
Fig.2.2)  was  derived  in  section  4.2.  It  is  expressed  in  Eq.(4.52)  as 


R2(x) 


(2L)2 

2k 


■p^x)- 

1  +a2 


(5.27) 


where 

F*!  (x) 

P,  (x)  =  -4-  (5.28) 

or 


and 

o 

a  =  .  (5.29) 

ai 

Since  sin_1(x)  is  continuous  and  single-valued  and  has  all  derivatives  on  an 
interval  including  x  =  0,  it  can  be  expanded  in  a  Mac  Laurin’s  series  at  x  =  0.  The 


47 


series  expansion  of  Eq.(5.27)  is 


R2(t) 


(2L)2 

2  K 


oo 


x  v 

k  =  o 


P1  (x) 
1  +a2 


^2k+i 


where  the  kth  coefficient  is 


1 


1  -3-5 . (2k  - 1 ) 

2-4-6 . (2k)  (2k+1) 


;  k  =  0 
;  k  >  o  . 


Substitution  of  Eq.(5.9)  into  (5.28)  results  in 


(5.30) 


(5.31) 


p!  CO 


j2rcfcx 

Re{R1  (x)  e  c  } 


j2JCf  x  .  .  -j27lf  x 
Ri  (x)  •  e  +  R-|  (x)  •  e 

2 


(5.32) 


Consequently,  the  series  expansion  for  R2(x)  becomes 


k  =  o 


1/(2  a2) 
1  +a2 


^2k  +  i 


j27lf  x  _  *  -j2jcf  x 

R-i  (x)  •  e  c  +  R-,  (x)  •  e  c 


2k+ 1 


(5.33) 


Following  the  same  development  used  to  obtain  Eq.(3.29)  from  Eqs  (3.11)  and 
(3.14),  it  follows  that 
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*2k+i 


FL(t)  = 


1\ 
k  =  0 


o  -,2k  + 1 

V(2o?)  ,2k 

- i-  -R,<t) 

1  +{X2 


j2Jtf  X 

Re{R1(x)  e  c  }+AR2(x)  (5.34) 


where  AR2(x)  denotes  the  terms  in  R2(x)  corresponding  to  the  spectral  zones 
other  than  those  centered  at  ±fc.  Define  the  complex  envelope 


p  oo  p2k+ 1 

A  ,  .  .  (2L)2  v  .  °k 

R*(t)L-—  A  k'lST 

c  k  =  0  d 


1  /  (2a2) 
1  +a2 


2k  +  i 


2k  ~ 

RlCO|  Rl(x)  •  (5.35) 


It  is  then  possible  to  express  the  autocorrelation  function  at  the  nonlinear  amplifier 
output  as 


1  -  i27tfcx 

R2(x)  =  ^Re{R2(x)  e  c}+AR2(x)  . 
d  Tc 


(5.36) 


The  bandpass  filter  following  the  nonlinear  amplifier  is  assumed  to  have  the  ideal 
transfer  function  whose  magnitude  is  shown  in  Fig. 2.4.  The  bandwidth  BRF  is  as¬ 
sumed  to  be  narrow  enough  such  that  only  the  spectral  zones  at  ±fc  are  passed 
with  negligible  distortion.  The  remaining  spectral  zones  are  assumed  to  be  per¬ 
fectly  attenuated.  It  follows  that  the  autocorrelation  function  at  the  RF  amplifier 
output  is  given  by 


1  -  i2Kicx 

R3(x)  =  ^Re {R2 (x) |f  e  c  } 

0 


(5.37) 
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5.3  Autocorrelation  Function  at  the  Transversal  Filter  Input 

From  Eq.(4.63)  the  complex  autocorrelation  function  at  the  transversal  filter 
input  is 


R2(x) 


R5(*)  = 


(5.38) 


Substitution  of  Eq.(5.35)  into  Eq.(5.38)  yields 


R5(t)  = 


o  oo  f'2k  +  1 

(2L)2  v  Ck 

2tc  ,2*  k  _2k  + 1 
k  =  o  2 


1/(20*) 
1  +  a2 


^2k+i 


|Ri  (x)|2k  Ri  (T) 


(5.39) 


It  should  be  recalled  that  the  above  result  is  valid  only  for  the  case  of  Gaussian 
signals  at  the  receiver  input. 


5.4  Autocorrelation  Function  Considerations  for  the  Ideal  Linear 
Receiver 


The  ideal  linear  receiver  provides  a  basis  of  comparison  for  performance  of 
the  nonlinear  receiver.  As  explained  in  Chapter  2,  the  receiver  nonlinear  behavior 
is  assumed  to  be  confined  to  the  RF  amplifier.  For  small  enough  inputs  the  nonlin¬ 
ear  amplifier  behaves  in  a  linear  fashion.  As  given  by  Eq.(2.3),  the  linear  gain  of 
the  RF  amplifier  is 


G  = 


(5.40) 


where  G  is  the  slope  of  the  nonlinear  input-output  characteristic  evaluated  at  the 
origin.  Consequently,  the  RF  amplifier  output  in  the  equivalent  ideal  linear  receiver 
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is  modeled  as 


j2Jtf  t 

r3  (t)  =  r2  (t)  =  G  •  r1  (t)  =  Re  {G  •  r1  (t)  •  e  c}  •  (5.41) 


It  follows  that  the  complex  envelopes  of  r3(t)  and  r2(t)  are  given  by 


r3  (t)  =  r2  (t)  =  G  •  ?!  (t)  .  (5.42) 

The  complex  autocorrelation  function  at  the  RF  amplifier  output  is 
R3(t)  =  R2  (x)  =  E  [?3  (t)  •  r3*  (t  — x)  ] 

=  G2  •  E  [r,  (t)  •  r/  (t  -  x)  ]  =  G2  •  R-i  (x)  .  (5.43) 

Using  the  result  of  Eq.(5.26),  we  obtain 

R3(t)  =  R2  (t) 


2G' 


j27CfnX 


20 


+  o_  e 


sin  (tiWt) 
7tWx 


.  (5.44) 


When  the  RF  amplifier  is  linear,  the  only  spectral  zones  in  Fig. 3.4  occur  at  ±fc. 
Hence, 


R2(x)  =R2(x)  . 

c 


(5.45) 


With  reference  to  Eqs  (4.63),  (5.44),  and  (5.45),  it  is  concluded  that  the  complex 
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auto-correlation  function  at  the  transversal  filter  input  is  given  by 


r5w 


R2(t) 

4 


G2 

2 


e 


j27IfDT 


20 


+  a  e 


w  2 
+  o 

n 


sin  (kWx) 

7tWx 


(5.46) 
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Chapter  6 


SIGNAL-TO-DISTURBANCE  RATIO  AT  THE 
TRANSVERSAL  FILTER  OUTPUT 

One  means  for  assessing  radar  performance  degradation  due  to  the  nonlin¬ 
ear  receiver  behavior  is  to  evaluate  the  signal-to-disturbance  ratio  at  the  transver¬ 
sal  filter  output.  Analytical  expressions  for  the  signal-to-disturbance  ratio  are 
obtained  in  this  chapter  for  both  the  linear  and  nonlinear  receivers  modeled  previ¬ 
ously.  The  transversal  filter  weights  are  then  optimized  to  maximize  the  signal-to- 
disturbance  ratios  for  the  following  cases  : 

Case  I  -  an  ideal  linear  receiver, 

Case  II  -  a  nonlinear  receiver  having  the  same  transversal  filter  as  was  used 
with  the  ideal  linear  receiver, 

Case  III  -  a  nonlinear  receiver  with  the  transversal  filter  optimized  to  account 
for  the  receiver  nonlinearity. 

Case  II  can  arise  when  the  design  engineer  either  does  not  know  of  the  nonlinear¬ 
ities  within  the  receiver  or  does  not  care  to  compensate  for  these  nonlinearities. 

6.1  Output  Stenal-to-Disturbance  Ratio  for  the  Ideal  Linear  Receiver 

Consider  the  ideal  linear  receiver  of  section  5.4.  A  block  diagram  of  the 
transversal  filter  is  shown  in  Fig. 2. 6.  Let  the  transpose  of  the  transversal  filter 
weight  vector  be 

WT  =  [w^g . wN]  b  (6.1) 

The  complex  envelope  of  the  transversal  filter  input  is  denoted  by  r5  (t) .  Let 

^  [  (N  -  k  +  1)Tg]  =  r^  ;  k  =  1,2 . N.  (6.2) 
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At  time  t  =  NTs,  the  complex  envelope  of  the  transversal  filter  output  is  then  given 
by 

r6(NTs)  =  (w  )T?5  (6-3) 

where 

^5  =  t  r51  r52 . •  (6.4) 

Noting  that  r6  (NTS)  is  a  scalar  quantity,  its  transpose  is  equal  to  itself.  Hence,  the 
average  power  in  r6  (NTg)  can  be  expressed  as 

E[fr6<NTs>|2]  “  Etf6(NTs>  '  V  <NTS)1 

-  e[?6(NTs)  ‘  <f6*  (NTs))T] 

=  e[(w  )Tf5- (r5")Tw] 

=  (w  )  E [ r5  -  (r5  )  Jw 
=  (w)TML-w  (6.5) 

where  the  correlation  matrix 

ml  =  e[f5-  (f5*)T]  (6.6) 

is  an  (N  x  N)  Hermitian  matrix.  The  superscript  L  is  used  to  indicate  that  the  anal¬ 
ysis  in  tnis  section  is  for  the  ideal  linear  receiver. 

The  entries  in  the  matrix  ML  are  readily  determined  from  the  complex  auto¬ 
correlation  function  r5  (t)  ,  as  is  shown  next.  By  definition, 
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E  (|V||  )  E  ^51  ^  ) 

E (^52^51  )  E(|^52|  ) 


E^5N^51  )  E^5N^52  ) 


The  uv,h  element  in  ML  is  given  by 


E(?51?5N  > 
E^52^5N  ) 


E(^5n|  ) 


(6.7) 


(ML)uv  =  E(r5„V>  =  E{"r5[(N-u  +  1)Ts]  r5*  [(N-v+1)T$]} 


=  R5[  (v-u)Ts]  .  (6.8) 

Recognizing  that 

R5[(u-v)Ts]  =  R5  [(v-u)TsJ  ,  (6.9) 

it  follows  that 

Rs(O)  r5(ts)  . r5[(N-dts] 

R5’  (Ts)  RS(0)  . RS[(N-2)TS] 


5  [(N-1)T  ]  r*  ( (N  — 2) T_] .  R5 (®) 

W 
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With  reference  to  Eq.(5.46),  note  that 


<ml), 


G2  r  2  j2rcfD(v-u)T 
=  R5[(v-u)Ts]  =  —  |^o2  e  D 


((v-u)Ts)' 


ac  6 


sin  (tcW(v-u)T  ) 

y2 - — 

n  tcW  (v  -  u)Te 


.(6.11) 


The  desired  signal  term  in  (ML)Uv  *s  readily  recognized  to  be 


2  2 

I  G  °s  j2rcfD(v-u)Tc 

<<)  =  ■  ©  D  s 


(6.12) 


whereas  the  disturbance  term  is 


((v-u)Ts)‘ 


<Md>, 


~2 

G  2 

—  a  e 
2  c 


2  sin  (tcW(v-u)Ts) 
*n  7tW(v-u)T_ 


(6.13) 


Consequently,  the  matrix  ML  can  be  decomposed  as 


ML  =  Mt  +  Mj 
s  d 


(6.14) 


where  (MLS  )uv  and  (MLd)uv  are  the  uv,h  elements  of  MLS  and  MLd,  respectively. 
When  Eq.(6.14)  is  substituted  into  Eq.(6.5),  the  average  power  in  r6(NT$)  be¬ 


comes 
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E[p6(NTs)|2]  =  (W  )TM^w+ (W  )TM^w  .  (6.15) 


The  signal  power  is  defined  to  be 

S  =  (w* )  TMg\iy  .  (6.16) 


Also,  by  definition,  the  disturbance  power  is  given  by 

D  =  (w  )TM^w  .  (6.17) 


As  a  result,  the  output  signal-to-disturbance  ratio  for  the  ideal  linear  receiver  is 

(w  )Vw 

(SDR)  j  =  - =— p—  •  (6.18) 

(W  )  Mjjir 


Having  determined  the  output  signal-to-disturbance  ratio  for  the  ideal  linear 
receiver,  the  nonlinear  receiver  is  treated  next.  Following  the  same  development 
as  in  section  6.1 ,  the  average  power  in  the  complex  envelope  of  the  transversal  fil¬ 
ter  output,  sampled  at  time  t  =  NTS,  is 


(6.19) 


where 


(6.20) 
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is  an  (N  x  N)  Hermitian  matrix.  The  superscript  NL  is  used  to  indicate  that  the 
analysis  in  this  section  is  for  the  nonlinear  receiver.  The  uvth  element  in  the  matrix 
MNLis 


(M  )uv  =  Rs[(v-u)TJ 


(6.21) 


where  R5  (x)  is  the  complex  autocorrelation  function  at  the  transversal  filter  input. 
From  Eq.(5.39) 


R5(t)  = 


n  00  p2k  +  1 

(2L)2  _ 

2jt  ,  2*  02k+i 

k  =  0  2 


1/(2o2) 
1  +a2 


2k+i 


—  gk  — 

-  Ri  (x)|  Ri(x) 


(6.22) 


In  terms  of  the  linear  gain  of  the  RF  amplifier,  note  that 

<2L)2  -  r2n2 

ST  =  G  °g  ■ 


(6.23) 


Hence,  R5  (x)  can  be  written  as 


00  Q2k  + 1 

R5(x)  =  G2<j2  y  - 

SL  «  k  p2k+1 
k  =  0  2 


1  /  (2a2) 
1  +a2 


2k  +  i 


-  2k  - 
|R-i  (X)|  •  Ri  (x)  . 


(6.24) 


From  Eq.(5.26) 


-  n  o  j27tfOT  2  2aw  2  sin  (»Wx) 

R,  (x)  =2  a2  •  e  +  cc '  8  "  +  <  ^7 — 


.  (6.25) 
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The  desired  signal  term  in  Rs  (t)  is  obtained  by  retaining  only  the  desired 
signal  term  in  R.|  (x) .  This  is  equivalent  to  evaluating  Rt  (x)  when  the  distur¬ 
bance  is  absent.  In  particular,  the  desired  signal  term  in  R1  (x)  is 


Ri(t) 


d(t)  =o 


2 

=  2cr  •  e 
s 


(6.26) 


Note  that 


-  2k  o  < 

Ri  (t)  .  4  A  =  (20 

d  (t)  =  o  s 


(6.27) 


It  follows  that  the  desired  signal  term  in  R5  (x)  is  given  by 


s2k+i 


«5(t)L,t>,n  =  QMlb 


■d  (t)  =  0  g  k  2k  +  i 
k  =  0  * 


1/(2ap 

1  +  a2 


2k+  1 


o  2k +1  j27CfDX 
<2°s>  •  e 


(6.28) 


As  a  result,  the  desired  signal  term  in  (MNL)UV  is 


(M?S 


00  p2k+ 1 

r2  9  V  h  .  J* _ 

g  ^  k  2k +  1 
k  =  o  * 


1/  (2a2) 
1  +a2 


2k +  1 


2  2k +  1  ]27ifD(v-u)Ts 


(2  V  • 6 


(6.29) 
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With  this  approach,  observe  that  (signal  x  signal)  cross  terms  which  arise  due  to 
the  nonlinearity  are  treated  as  desired  signal  components  since  they  all  contribute 
to  the  Doppler  component  evaluated  at  t  =  (v-u)Ts. 

As  a  check,  it  is  now  shown  that  the  k  =  0  term  reduces  to  Eq.(6.12)  for  the 
linear  receiver.  From  Eqs.(3.16)  and  (5.31)  note  that 

Cj  =  b0  =  1  .  (6.30) 


Hence,  the  k  =  0  term  in  Eq.(6.29)  is 


(MNL) 


UV 


gV 

1/  (2o2) 

k.o"  2 

1  +  a2 

]27CfD(v-u)T 

2a“ • e  u  s 
s 


(6.31) 


However,  a  is  defined  in  Eq.(4.12)  as  the  ratio  of  cg  to  a-,.  Therefore, 


<L> 


UV 


gV 


g.  1  .  qi  ;o2  J2*'p <v-“) ts 

S 


k  =  0  2  2a2  o2  +  a2 


.  (6-32) 


The  receiver  behaves  in  a  linear  fashion  when  ai2«  og2.  Using  this  inequality,  it 
follows  that 


(Msnl> 


UV 


k  =  0 


gVs  j2^D(v-u>TS  fML. 
—  9  =(MS)uv' 


(6.33) 


As  expected,  the  linear  result  is  a  special  case  of  the  nonlinear  result. 

The  disturbance  term  in  (MdNL)uv  is  simply  obtained  by  substracting 
(Msnl)uv  from  (Mnl)uv  .  This  is  valid  because  (MNL)UV  is  composed  of  (signal  x 
signal)  cross  terms  plus  (signal  x  disturbance)  cross  terms  plus  (disturbance  x  dis¬ 
turbance)  cross  terms.  By  subtracting  the  (signal  x  signal)  cross  terms, we  are  left 
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with  the  disturbance  terms  which  are  given  by 


(MdL>„u=  <MNL>uv-<MsL>  • 


UV 


(6.34) 


uv 


Let  (MSNL)UV  and  (M<jNL)uv  be  the  uvth  elements  of  the  matrices  MSNL  and 


MdNL,  respectively.  Then 


MNL=mNL+MNL 
s  d 


(6.35) 


Substitution  of  Eq.(6.35)  into  Eq.(6.19)  results  in 

E[|r6(NTs)|2J  =  (w*)TMg  Lw  +  (w  )TM^Lw 


The  signal  and  disturbance  powers  are  defined  to  be 


and 


(6.36) 


S  =  (w*)TMg  LW 

(6.37) 

D  =  (w  )TMdNLw  . 

(6.38) 

Consequently,  the  output  signal-to-disturbance  ratio  for  the  nonlinear  receiver  is 


(SDR)nl  = 


(W  ^ 

(W  )TM^LW 


(6.39) 


The  uvth  elements  in  the  matrices  MNL,MSNL  and  MdNL  consist  of  infinite  se¬ 
ries.  Evaluation  of  these  elements  raises  the  question  of  how  many  terms  should 
be  summed  in  order  to  obtain  an  adequate  approximation.  From  Eqs  (6.21)  and 
(6.24) 
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2k +  i 

•  RiI(v-u)Ts]  .  (6.41) 

In  addition,  with  reference  to  Eq.(5.26), 

Ri  [(v-u)Ts]|<2[a|  +  a2  +  o2]  .  (6.42) 

Because  the  desired  signal,  clutter  and  noise  are  assumed  to  have  zero  means, 
their  variances  are  given  by 

RS(0)  =  as’Rc(0)  =  ac,Rn(0)  =  an  -  (6.43) 

From  Eq.(5.4)  it  follows  that  the  variance  of  the  input  signal  rj(t)  is 

°i  =  R-i  (0)  =  °s  +  °c +  °n  *  (6.44) 
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Consequently,  the  inequality  in  Eq.(6.42)  becomes 


Ri  [  (v-u)T  ]  <2a* 


(6.45) 


Use  of  Eq.(6.45)  in  Eq.(6.41)  results  in 


.2k+i 


/  fcyjN  L\  2  Y  h  k  [” 2 

(  uv  9k?o  k'  22kt1  il+a2. 


2k+i 


(6.46) 


Relative  to  the  infinite  series  in  Eq.(6.46),  define  the  finite  series 


>2k+i 


2k+i 


(6.47) 


Values  of  T(M)  are  tabulated  in  Table  6.1  as  a  function  of  a  and  M.  a  =  0.05  corre¬ 
sponds  to  the  situation  where  the  rms  value  of  the  input  equals  twenty  times  the 
linear  range  of  the  nonlinearity,  as  determined  in  Eq.(2.7).  Assuming  it  is  unlikely 
that  the  nonlinearity  will  be  driven  harder,  we  only  consider 


a  2  0.05  ■ 


(6.48) 


Because  the  kth  term  in  T(M)  decreases  with  increasing  values  for  a,  the  rate  of 
convergence  of  T(M)  is  most  severe  for  a  =  0.05.  When  a  =  0.05,  it  is  seen  from 
Table  6.1  that  T(M)  has  almost  converged  for  M  =  11.  The  relative  change  that  re¬ 
sults  when  the  next  term  is  added  is 


T (12)  -T (11) 
T  ( 1 1 ) 


=  I.SxlO-'3  • 


(6.49) 


Ct  =  0.05 
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M 

a  =  0 

a  =  0.05 

a  =  0.1 

ct  —  0.5 

a=  1 

0 

0.50000 

0.498753 

0.495049 

0.400000 

0250000 

1 

0.56250 

0.560786 

0.555710 

0.432000 

0.257810 

2 

0.58590 

0.583896 

051191  A 

0.439667 

0258541 

3 

0.59810 

0,595884 

0 .589353 

0.442225 

0258636 

.4 :  ■ 

0.60557 

0.603188 

0.596184 

0.443227 

0.258650 

5 

0.61061 

0.608098 

0.600706 

0.443660 

0.258653 

6 

0.61425 

0.611616 

0.60390 

0.443860 

0.258653 

7 

0.61699 

0.614258 

0.606262 

0.443956 

0.258653 

8 

0.61913 

0.616311 

0.608071 

0.444005 

! 

9 

0.62085 

0.617951 

0.609495 

0.444029 

| 

10 

0.62226 

0.619290 

0.610640 

liliii 

11 

0.62344 

0.620403 

0.611577 

0.444049 

12 

0.62444 

0.621342 

0.612356 

0.444053 

! 

13 

0.62530 

0.622144 

0.613012 

0.444055 

I 

14 

0.62604 

0.622836 

0.613570 

0.444056 

| 

15 

0.62669 

0.623440 

0.614049 

0.444057 

i 

16 

0.62727 

0.623970 

0.614464 

0.444057 

i 

Table  6. 1  -  Values  of  T(M)  as  a  Function  of  a  and  M 

A  plot  of  T(M)  as  a  function  of  M  for  values  of  a  ranging  from  0.05  to  1  is  shown  in 
Fig-6.1.  It  is  concluded  that  1 2  terms  are  adequate  for  evaluating  the  elements  in 
Mnl.  It  is  readily  shown  that  this  conclusion  also  holds  for  calculating  the  terms  in 
Msnl  and  MdNL. 
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Figure  6. 1  -  Plot  of  T(M)  as  a  function  of  M  for  values  of  a 

ranging  from  0.05  to  1. 


6.3  Determination  of  the  Transversal  Filter  Weights 

The  weights  of  the  transversal  filter  are  selected  so  as  to  maximize  the  sig- 
nal-to-disturbance  ratio  at  the  transversal  filter  output.  As  mentioned  previously, 
there  are  three  cases  to  consider.  Each  is  discussed  separately. 


Case  I 

This  case  deals  with  the  ideal  linear  receiver.  From  Eqs  (6.18)  and  (6.12)  the 
output  signal-to-disturbance  ratio  is 


(SDR) 


-H.j- 

yy  MgW 

1  =  - H.  J  - 
W  Mjjyy 


(6.50) 


where  the  superscript  H  denotes  conjugate  transpose  and  the  uvth  element  of  MSL 
is 
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(6.51) 


L  G  °s  j2*fD(v-u)T 
(Mb  =— D  s 


For  convenience,  let 


ak  =  e 


j2*fDkTs 


;  k  -  1 ,2 . N  . 


(6.52) 


Define  the  vector  a  such  that 


3  =  [3^2 . aN] 


(6.53) 


Note  that  the  matrix  MSL  can  be  expressed  as 


r2  2 

ML  s  a*  J 

=  — - —  •  a  a 

s  2 


(6.54) 


If  follows  that 


G2a2  -H  *  T-  G2a2  -H  *|2 
(SDR)j-  •  _H  -  '  -HmL- 


(6.55) 


Let  x  and  y  denote  two  N-component  complex  vectors.  Because  Md  L  is  a 
positive  definite  (N  x  N)  Hermitian  matrix,  it  can  be  shown  that 


(s.y)  =  *HM[jy 


(6.56) 


is  an  inner  product  of  the  vector  space  of  N-component  complex  vectors.  An  im¬ 
portant  property  of  the  inner  product  is  the  Schwartz  inequality 


(x,y)|2^  (x.x)  (y,y) 


(6.57) 
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where  equality  holds  if  and  only  if 


i  =  ky 


(6.58) 


and  k  is  an  arbitrary  complex  constant  [2].  In  terms  of  the  inner  product,  note  that 


WHS*  =WHM^(M^)  a*  =  (W,  (M|p  a*)  .  (6.59) 


Also, 


H  I 

W  MJjw  =  (w,w) 


(6.60) 


As  a  result,  the  output  signal-to-disturbance  ratio  can  be  written  as 


gV 


(SDR)  j  = 


(w.  (M|p  a* ) 


~H.  ,L  - 
W 


(6.61) 


Application  of  the  Schwartz  inequality  yields) 


G2<j?  (w,w)  ((mJ-)  S,  (M*5>  8) 
(SDR),*;-^ - 0  0 


(w,  w) 


G2o2 


^  s  i  -1  .  i  -i  . 

£  — 2~  •  ( (Mjp  a,(Mjp  a) 


(6.62) 


Equality  holds  when 


-  i  -i  * 

w  =  k-  (M^)  a  . 


(6.63) 
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Because 


-1  _1 
((M|j)  >  =  <Md>  • 


(6.64) 


the  maximum  signal-to-disturbance  ratio,  which  is  obtained  when  the  weight  vec¬ 
tor  is  chosen  according  to  Eq.(6.63),  is  given  by 


gV 


(SDR)  j  =  — ^  •  ((M|j)  s',  (M|j)  s') 


G2o2  H 

-r?-8T((M^)  )  M{j(M{j)  s' 


G2<j? 


2J8T<mL)'  s’ 


(6.65) 


To  verify  that  this  is,  indeed,  the  signal-to-disturbance  ratio  when  Eq.(6.63)  is  val¬ 
id,  the  weight  vector  is  substituted  into  Eq.(6.55).This  yields 


-1  n  -1 

gV  k"  3T((Mj)  )  aVk(Mt)  a* 

(SDR),  =  V - ^-TT-H - V- 

u*  «T / iL\  \  kjLiwkjL\  * 


((M|p  )  MJjk(MJp  a 


G2o2 


S  T  |  —  1  * 
~2~ '  aT(Mjj)  a 


(6.66) 


The  result  obtained  in  Eq.(6.66)  is  consistent  with  Eq.(6.65). 
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Case  II 


In  this  case  the  receiver  is  assumed  to  be  nonlinear.  However,  the  presence 
of  the  nonlinearity  is  unknown  to  the  system  designer.  Hence,  the  “optimum" 
weight  vector  is  chosen  as  though  the  receiver  is  linear.  With  reference  to 
Eq.(6.39),  the  signal-to-disturbance  ratio  is  given  by 


whmsnlw 

<SDR,n  ■  Th^UT 


(6.67) 


where  the  weight  vector  is  determined  according  to  Eq.(6.63).  This  weight  vector 
does  not  result  in  a  maximum  signal-to-disturbance  ratio  because  it  depends  upon 
the  disturbance  covariance  matrix  of  the  linear  receiver  instead  of  the  disturbance 
covariance  matrix  of  the  nonlinear  receiver. 


Case  III 

Once  again,  the  receiver  is  assumed  to  be  nonlinear.  The  weight  vector  for 
this  case  is  chosen  to  maximize  the  signal-to-disturbance  ratio.  From  Eqs  (6.39) 
and  (6.29)  the  output  signal-to-disturbance  ratio  is 


(SDR) 


WHMg  Lw 

w^M^w 


(6.68) 


where  the  uvth  element  of  MJ-  is 


<^L> 

»  uv 


-  ,^2 


>2k  +  i 


G  cr  Y  bk  • 

g k  2k +  i 
k  =  0  2 


1/(2op 

1  +  a2 


2k  + 1 


•  (2<5> 


2k*1  T  (6^69) 


69 


Let  A  denote  the  constant 


»2k+ 1 


>2  2  x~<  i  k 


A  =  GV  X  b 


g,^  k  2k  +  i 
k  =  o  2 


1/(2  o‘) 
1  +  a2 


^2k  +  i 


(2o|) 


2k+i 


(6.70) 


In  terms  of  the  vector  a  defined  by  Eqs  (6.52)  and  (6.53)  the  matrix  MSNL  can  be 
expressed  as 


=  A-3*3T  . 


Hence,  the  signal-to-disturbance  ratio  becomes 


-H  *  T-  -H  * 

fcnn>  _aW33  W_.  W  3 

(SDR)I»_A  "  A 


(6.71) 


(6.72) 


w  MjJSy 


w  M^w 


Because  M<jNL  is  a  positive  definite  Hermitian  matrix,  it  is  possible  to  define  the  in¬ 
ner  product  by 


(x,y)  =  xHM^Ly  . 


(6.73) 


Following  the  same  procedure  as  in  Case  I,  it  can  be  shown  that  the  optimal 
weight  vector  is  given  by 

W  =  k-(Mjpa  (6  74) 

where  k  is  an  arbitrary  complex  constant.  The  maximum  signal-to-disturbance  ra¬ 
tio  is 

(SDR) in  =  A-3T(m|p  V  .  (6.75) 
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The  analytical  results  provided  in  this  chapter  are  used  in  the  following 
chapters  to  calculate  the  signal-to-disturbance  ratios  for  the  three  cases  consid¬ 
ered. 
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Chapter  7 

FREQUENCY  DOMAIN  ANALYSIS 


The  previous  chapters  utilized  a  time-domain  approach  to  develop  the  sig- 
nal-to-disturbance  ratio  at  the  transversal  filter  output.  However,  much  physical  in¬ 
sight  can  be  gained  by  interpreting  the  results  in  the  frequency  domain.  Assuming 
that  the  input  signals  are  sampled  in  time,  the  continous  Fourier  transform  of  a 
discrete-time  series  is  first  discussed.  This  leads  to  development  of  the  transver¬ 
sal  filter  transfer  function.  With  respect  to  the  power  spectral  densities  of  the  clut¬ 
ter,  thermal  noise,  and  transmitted  signal,  the  effect  of  processing  a  finite  number 
of  samples  is  considered.  Finally,  a  relation  is  derived  for  the  power  spectral  den¬ 
sity  at  the  transversal  filter  output. 

7.1  Continuous  Fourier  Transform  of  Discrete  Time  Series 

Let  x  [nT  ]  denote  a  complex  valued  discrete  time  sequence.  Its  continu- 

5 

ous  Fourier  transform  is  defined  [8]  to  be 

00  -j2itnfT 

X  (f)  =  £  x[nTs]e  s.  (71) 

n  =  -oo 


Note  that 


~  -j2*n(f±Ajrs 

1  *(nTs]e 


=  X  *lnTS)0 

n  =  -oo 


-]27tnfT 


seTj2itnk  =  x  (()  .  (7.2) 

r 


Hence,  X  (f)  is  periodic  with  period  1/TS.  The  inverse  Fourier  transform  is  given 
by 
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X 


[nTs]  =  Ts 


]27tnfT 
e  sdf  . 


(7.3) 


It  is  convenient  to  suppress  dependence  on  the  sampling  interval  Ts  by  introduc¬ 
ing  the  normalized  frequency  variable 


Q  =  2n 


-  2*,Ts  • 


(7.4) 


Note  that  £2  is  in  radians/second  whereas  f  is  in  Hertz.  In  terms  of  £2,  Eqs  (7.1) 
and  (7.3)  become 

oo 

X(Q)  =  £  5[nTs]e“,nQ  (7.5) 

n  =  -oo 


and 

71 

x  [nTsl  =  ±  J  X  (Q)  einnd«  .  (7.6) 

-71 


Observe  that 

X(Q±27ck)  =  ]T  x[nTs]e^n(Q±27lk) 
n  =  -oo 

=  2  i(nTs]e"inVi2,ink  =  X(£J)  .  (7.7) 


Consequently,  X  (Q)  is  periodic  with  period  2tt. 
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7.2  TYansfer  Function  of  the  Transversal  Filter 


As  given  by  Eq.(2.9),  the  transversal  filter  output  is  determined  by 

N 

r6(t)  =  I  wk#?5[t-(k-1)Ts] 

k  =  i 


(7-8) 


where  r5  (t)  denotes  the  transversal  filter  input.  Assuming  that  the  signals  are  uni¬ 
formly  sampled  with  sampling  interval  Ts,  the  filter  output  at  the  nth  sampling  in¬ 
stant  is 

N 

r6[nTsl  =  Z  wk* f5[(n-k+1)Ts]  .  (7.9) 

k  =  i 


Let  X5  (£2)  and  X6  (£2)  denote  the  continuous  Fourier  transforms  of  r5  (nTg)  and 
r6(nTs) ,  respectively.  Using  Eq.(7.9)  in  Eq.(7.5)  and  interchanging  the  order  of 
summation  results  in 


x6(n>  =  Z  tnTsi e  Jn^  =  Z 


n  =  -oo 


N 

I  "kl 

k  =  1 


n  =  -oo 


N 

X  wk‘r5[(n-k+1)TsJ  [e 

Ik- 1 


.-jn£2 


oo 

X  r5((n-k+1)T  ]e~jnQ 


n  =  -oo 


(7.10) 


Introduce  the  change  of  variables 


I  =  n-k+ 1  . 


(7.11) 


It  follows  that 


X6(Q)  =  X  *kVi(k“1)n 

k  =  i 


I  ^5 1*1.1  e 


-jin 


I  =  -OO 
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(7.12) 


N 

=  x5(Q)  2 

k  =  i 


w„  e 


-j(k-  i)Q 


By  definition,  the  transversal  filter  transfer  function  is 


Hjp(Q) 


^6 (^)  -  *  -j(k-i)Q 

^ -  =  X  wk  e 

Xs(O)  k=  i 


(7.13) 


As  explained  in  Chapter  6,  the  weights  wk  are  determined  so  as  to  maximize  the 
signal-to-disturbance  ratio  at  the  filter  output. 


7.3  Power  Spectral  Density  of  Complex  Envelope  at  Receiver  Input 

Assuming  wide-sense  stationary,  consider  the  continuous  autocorrelation 
function  R  (t)  where  x  denotes  the  difference  between  two  time  instants.  For  uni¬ 
form  sampling,  with  a  sampling  interval  Ts,  the  time  difference  is  an  integer  multi¬ 
ple  of  Ts.  Consequently,  we  let 

t  =  kT$  (7.14) 

where  k  is  an  integer.  The  associated  power  spectral  density  S  (Q)  is  the  contin¬ 
uous  Fourier  transform  of  the  sequence  R  (kT  ) .  With  reference  to  Eqs.(7.5)  and 

w 

(7.6),  this  results  in  the  Fourier  transform  pair 

OP 

S(Q)  =  £  R[kTs]e'jkQ  (7.15) 

k  =  -oo 


and 

7t 

R[kTs]  =  jL  |S(fl)9ikQdn  .  (7.16) 

-n 

Given  that  a  target  is  present,  the  received  signal  consists  of  the  sum  of  the 
signal  reflected  from  the  target ,  clutter  and  thermal  noise.  Because  the  transver- 
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sal  filter  processes  N  samples  at  a  time,  the  target,  clutter,  and  thermal  noise  se¬ 
quences  are  effectively  truncated  to  N  samples.  In  this  section,  we  study  the  effect 
of  this  truncation  on  the  corresponding  power  spectral  densities. 


7.3.1  Effect  of  Truncation  on  the  Clutter  Power  Spectral  Denitv 

The  autocorrelation  function  of  the  clutter  complex  envelope  is  given  by 
Eq.(5.22)  as 


Rc(x>  =  2ac'e 


(7.17) 


For  the  sampled  clutter,  the  corresponding  correlation  sequence  is 


k2n 


Rc[kT  ]  =  2<j2  -e 


20 


w 


(7.18) 


For  convenience,  define  the  clutter  correlation  parameter  pc  to  be 


RC[TS] 

Pr  =  - - 

c  FUO] 


2o2  •  e 


20 


w 


20 


=  e 


w 


2  o! 


(7.19) 


Observe  that 


0£pc£1 


(7.20) 


In  terms  of  pc  and  Rc  [0] ,  the  clutter  correlation  sequence  becomes 


Rc[kT  ]  =  Rc[0]  •  p*  . 


(7.21) 
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The  clutter  samples  are  said  to  be  highly  correlated  when  pc=l  and  are  said  to  be 
uncorrelated  when  pc=0. 

A  sketch  of  the  complete  correlation  sequence  is  shown  in  Fig.7.1.  When  N 
clutter  samples  are  processed,  the  maximum  time  displacement  between  the 


samples  is  (N-l)Ts.  Hence,  the  allowable  values  for  k  are  k=0,±1,±2, . ±(N-1). 

The  effect  is  to  truncate  the  correlation  sequence,  as  is  shown  in  Fig. 7.2.  The 
truncated  correlation  sequence  is  defined  to  be 


{Rc[kTg]  ;  -(N-1)£k<(N-1) 

(7.22) 

0  ;  otherwise  . 

The  power  spectral  density  corresponding  to  the  truncated  correlation  sequence 
is  given  by 

(N-i)  ko 

ScT(Q)  =  I  Rc[kTsJe  ■  (7.23) 

k  =  -(N-i) 
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-<N-1)  0  <N-1) 

Fiaure7.2  -  Sketch  of  the  truncated  correlation  sequence  .RcitkTs^ 

The  truncated  correlation  sequence  may  no  longer  be  a  valid  correlation  se¬ 
quence  in  the  sense  that  the  corresponding  power  spectral  density  may  assume 
negative  values.  This  is  illustrated  in  Fig. 7.3  where  ScT  (Q)  is  plotted  for  various 
values  of  pc  and  N.  For  convenience,  the  frequency  axis  is  normalized  by  tc  so 
that  the  variation  in  Q  from  -  n  to  7t  is  transformed  into  a  variation  of  QJn  from  -1 
to  1 .  When  pc=0 


RdlkTsl- 


{ 


2c" 


;k=0 
;k*0  . 


(7.24) 


Hence, 


Sct(Q)=2<j2  (725) 

independent  of  the  choice  of  N.  This  constant  power  spectral  density  is  shown  in 
Fig.  7.3.(a).  For  pc=0.8  and  N=5,  ScT(Q)  is  essentially  Gaussian  shaped  and 
nonnegative,  as  is  shown  in  Fig.7.3.(b).  However  for  pc=0.9  and  N=5,  negative 
side  lobes  are  observed  in  SCT(Q).  This  is  shown  in  Fig.7.3.(c).  As  seen  in 
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Fig.7.3.(d),  the  negative  side  lobes  disappear  when  N  is  increased  to  10  and 
pc*0.9. 


Figure  7.3  -  Plots  of  SrT  (Q)  vsQ/n. 


(a)  0^=0  and  N=5.  (b)  oz=0.8  and  N=5,  (c)  oz=0.9  and  N=5, 

(d)  o^=0. 9  and  N=  1 0. 


Because  we  are  interested  in  those  cases  for  which  the  clutter  is  much  larg¬ 
er  than  both  the  signal  reflected  from  the  target  and  the  additive  thermal  noise,  the 
clutter  spectrum  dominates  and  N  is  chosen  large  enough  so  as  to  avoid  negative 
side  lobes  in  ScT  (Q) .  By  choosing  N  large,  the  edge  effects  of  the  truncated  cor¬ 
relation  sequence  are  minimized.  Define  the  ratio  of  the  smallest  to  largest  values 
of  RCT(kTs^  t0  be 


UN  = 


RctKN-DTsI 

RcTl°l 


2o2  •  p(N-1)‘ 

c  c _ 

2o? 


=  o(N-1)‘ 


(7.26) 
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Solution  for  N  results  in 


N  =  1  + 


iog(UN) 

a|  log  (pc) 


(7.27) 


In  order  to  approximate  Rc  (kTg)  as  closely  as  possible,  UN  is  chosen  equal  to 
10‘9.  The  values  of  N  needed  for  Un=  10"9  are  plotted  in  Figures  7.4  through  7.6 
for  various  values  of  pc.  In  Figure  7.4,  0<pc<1.  Because  of  the  large  range  of  N, 
expanded  plots  are  shown  for  0<pc<0.9  in  Fig. 7. 5  and  for  0.9<pc<1  in  Fig.7.6.  In 
our  analysis  the  maximum  values  allowed  for  pc  is  pc=0.9.  From  Fig. 7.5,  the  re¬ 
quired  value  of  N  is  N=15.  As  a  result,  it  is  assumed  in  our  analysis  that  the  trans¬ 
versal  filter  processes  no  fewer  than  1 5  samples. 


Figure  7.4  -  Values  of  N  needed  for  Um=  1 0^.  0<o^  1. 
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7.3.2  Effect  of  Truncation  on  the  Thermal  Noise  Power  Spectral  Density 

From  Eq.(5.25),  the  power  spectral  density  of  the  thermal  noise  complex  en¬ 
velope  is 


RnW 


2  S'n  (icWt) 
°n  TtWx 


(7.28) 


where  W  is  the  bandpass  bandwidth.  Assuming  uniform  sampling  of  the  noise, 
with  the  sampling  interval  Ts,  the  corresponding  correlation  sequence  is 


Rn[kTs]  =2o2n- 


sin  (7tWkTs) 
rcWkT 


(7.29) 


Because  the  magnitude  of  sin(7tWkTs)  is  bounded  by  unity,  the  bound  on  the  mag¬ 
nitude  of  Rn  (kTs)  for  k?to  is 


Rn  (kTsJ 


(7.30) 


Recall  that  Ts  is  equal  to  the  time  interval  between  consecutive  pulses  of  the 
transmitted  pulse  train.  Since  this  is  typically  on  the  order  of  milliseconds  while  W 
is  on  the  order  of  Megahertz,  it  is  assumed  that 


rcWkTs|»2a2  ;k  =  ±l,±2, .  (7.31) 


Thus,  Rn  (kTg)  is  approximated  by 


Rn[kTs] 


0 


;k*0 


(7.32) 
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A  sketch  of  the  approximate  autocorrelation  sequence  of  the  sampled  ther- 
mal  noise  is  shown  in  Figure  7.7.  Because  of  this  approximation,  the  power  spec¬ 
tral  density  of  the  truncated  noise  sequence  is 

(N--I) 

SnT(Q)  =  X  Rn[kTs]e-JkSi-2^  (7.33) 

k  =  -(N-i) 


Fiaure7. 7  -  Approximate  autocorrelation  sequence  of  the  sampled 

thermal  noise. 

independent  of  the  choice  of  N.  It  is  concluded  that  truncation  has  a  negligible  ef¬ 
fect  on  the  sampled  noise  spectrum. 

7.3.3  Effect  of  TVuncation  on  the  Desired  Siena!  Power  Spectral  Density 

The  autocorrelation  function  of  the  desired  signal  complex  envelope  is  given 
by  Eq.(5.16)  as 

_  j2ftfnT 

Rs(x)=2o|e  D.  (7.34) 

For  the  sampled  signal,  the  associated  correlation  sequence  is 

_  j27tfnkT 

Rs[kTs)  =2o|.e  D  .s  (7.35) 
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With  reference  to  Eq.(7.15),  the  power  spectral  density  of  the  desired  signal  is 


Ss(Q)  =  £  Rs  fkTs] e_)k£1  =  2<5s '  1  8 

k  =  -oo  k  =  — 


i2*fDkTs  -jkQ 

®  .  (7-36) 


For  convenience,  let 


QD  =  2lt,DTs 


(7.37) 


The  expression  for  the  power  spectral  density  then  simplifies  to 

-jk(Q-QD) 


Ss(fl)  =  2o|-  X  e 

k  =  -<= 


(7.38) 


Further  simplification  is  possible  by  using  the  Poisson  sum  formula  [6] 

00  -jk(a-Qn) 

£  e  u  =  2tc  £  8(Q-QD-2nji)  (7.39) 

k  =  oo  k  =  -oo 

where  5(.)  denotes  the  usual  Dirac  delta  function.  Use  of  Eq.(7.39)  in  Eq.(7.38)  re¬ 
sults  in 


oo 

Ss(Q)  =  47ia|  £  8(C2-QD-2n7i)  .  (7.40) 

k  =  -oo 

The  power  spectral  density  is  seen  to  consist  of  a  periodic  impulse  train,  as  shown 
in  Fig. 7.8. 

The  power  spectral  density  of  the  truncated  signal  sequence  is 

(N-i)  -jk<a-oD) 

SsJ(Q)=2c2s-  £  e  .  (7.41) 

k  =  -(N-i) 
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Figure  7.8  -  Power  spectral  density  of  desired  signal. 

This  is  recognized  as  the  sum  of  a  geometric  progression  which  equals 


S8T(G> 


=  2  a 


-j(fl-aD)  -j(N-i)<n-nD)  -j(N-i) 

2  e _ e _ -e _ 

-j(Q-Qn) 
e  u  -1 


-J 


Q-Q 


u  -i(N-l)(n-nn)  j(N-l)(n-Qn)^ 


e 


e 


D  •*'  "  2' 

-  e 


2o! 


V 


i(^) 

e  e  -e 


,Q-Qd^  YQ~Qnv\ 


(7.42) 


The  zeros  of  the  numerator  occur  when 


Consequently,  the  maxima  of  Ssj(Q)  occur  for 


(7.46) 


(7.47) 


The  power  spectral  density  of  the  truncated  signal  sequence  is  sketched  in 
Fig.7.9.  Due  to  truncation,  negative  side  lobes  are  seen  to  exist  about  the  main 
lobes.  These  can  be  eliminated  by  windowing  the  data.  However,  this  is  not  in¬ 
cluded  here  because  the  time-domain  analysis  of  the  previous  chapters  does  not 
involve  windowing.  In  addition,  even  with  negative  side  lobes,  the  area  under  the 
signal  power  spectral  density  equals  the  average  signal  power  as  evaluated  in  the 
time  domain. 
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Neverthless,  it  is  instmctive  to  examine  how  the  choice  of  N  affects  the 
shape  of  the  signal  power  spectral  density.  With  reference  to  Eq.(7.42),  side  lobe 
maxima  closely  coincide  with  peaks  in  the  numerator  due  to  the  relatively  flat 
shape  of  the  denominator  as  a  function  of  Cl  Peaks  of  the  numerator  result  when 

^^(Q-n0)  =  (2k+i>| 

or 

(fl-Q  )  =  ;k -0+1+2 .  (7.48) 

'  D'  2N-1 

For  this  value  of  (Q-ftD),  argument  of  the  sine  function  in  the  denominator  is 


Ssy  (Q) 


Figure  7.9  -  Power  spectral  density  of  the  truncated  signal  sequence. 


Im-Q  )  =  2k  +  1  71 

2  2N-1 


Letting  k=1 ,  the  amplitude  of  the  first  side  lobe  is 


(7.49) 


2  a* 


A,  = 


.  ,3  *  , 

S,n  2  2N^-T 


(7.50) 
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The  ratio  of  A1  to  the  amplitude  of  the  main  lobe  is 


r 


Ai 

2c*(2N-1) 


1 


3  K 

(2N-1) 

Sin<22N-t) 

(7.51) 


A  plot  of  r  as  a  function  of  N  is  shown  in  Fig.  7.10.  For  N  >5,  note  that  r  approach¬ 
es  its  asymptotic  value  of 


2  (2N  -  1 ) 
(2N- 1)3ti 


=  0.212 
3ji 


(7.52) 


Figure  7.10-  Plots,  as  a  function  of  N.  of  1)  Ratio  of  first  side  lobe 
amplitude  to  that  of  the  main  lobe  and  2)  width  of  the  main  lobe. 


The  width  of  the  main  lobe  at  QD  is  defined  to  be  the  frequency  spacing  be¬ 
tween  the  two  zeros  of  the  power  spectral  density  located  at 


Q 


1 


2jc 

2N^i 


and 


=  QD  + 


2k 

2N-  1 
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Thus,  the  main  lobe  width  is  given  by 


4  k 

2N-  1 


(7.53) 


Fig.7.10  shows  a  plot  of 


2k 


2 

2N-1 


(7-54) 


as  a  function  of  N.  Observe  that  the  width  approaches  zero  as  N  becomes  large. 
For  N=1 5, 

Q  -Q 

— 1  =  0.069  .  (7.55) 


The  plots  in  Fig.7.10  are  expanded  in  Fig. 7. 11  for  1<N<20. 


Figure  7.11  -  Expanded  plots,  as  a  function  of  N.  of  1)  Ratio  of  first 
side  lobe  amplitude  to  that  of  the  main  lobe  and 
2)  width  of  the  main  lobe. 

Plots  of  the  signal  power  spectral  density  for  N=5,1 0,1 5,50  and  Qd=0  are 
presented  in  Figures  7.12  *  7.15,  respectively. 


Magnitude 


Fiaure  7.12-  Discrete  sianal  cower  spectral  den  si 


for  N=5  andQn=0. 


Fiaure  7.13-  Discrete  sianal  power  spectral  densit 


for  N=1 0  and  Qn=0. 


Magnitude 


Figure  7. 14  -  Discrete  signal  power  spectral  density 
for  N=1 5  and  Qn=0. 


for  N=50  and  Qn=0. 
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7.4  Power  Spectral  Densities  at  the  Output  of  the  Nonlinearity 

As  explained  in  section  6.2,  each  element  of  the  correlation  matrix  of  the  sig¬ 
nal  at  the  nonlinearity  output  consists  of  (signal  x  signal)  cross  terms,  (signal  x  dis¬ 
turbance)  cross  terms  and  (disturbance  x  disturbance)  cross  terms.  By  separating 
out  the  (signal  x  signal)  cross  terms,  it  is  possible  to  express  the  correlation  matrix 


MNL  =  mN1-  +  mNL 
s  d 


(7.56) 


where  L  is  the  signal  correlation  matrix.  The  disturbance  correlation  matrix, 
m[Jl  ,  contains  the  (signal  x  disturbance)  cross  terms  plus  the  (disturbance  x  dis¬ 
turbance)  cross  terms  and  is  obtained  by  substracting  m^l  from  mnl. 

5 

From  Eq.(6.29),  the  uvth  element  of  L  is  given  by 


<msL) 

UV 


=  5V0IV 


*2k+i 


,2k  + 1 


1  +a2 


2k +  i 


9  2k  +  i  j2rcfn(v-u)T 
(2o‘)  e  u  s 


(7.57) 


This  element  corresponds  to  a  lag  of  (v-u)Ts  seconds.  As  a  result,  the  power  spec¬ 
tral  density  of  the  desired  signal  at  the  output  of  the  nonlinearity  can  be  obtained 
from 


§sT<«)  =  S  <MsNL)  e-i(,-u,£i+  £  (m"l) 

b  U1  .  5  tv 

u  =  N  v  =  i 


(7.58) 


The  summation  over  u  involves  elements  from  the  first  column  of  ,  incorporat¬ 
ing  the  lags  extending  from  -(N-1  )TS  to  -Ts,  while  the  summation  over  v  involves 
elements  from  the  first  row  of  M'  ,  incorporating  the  lags  from  0  to  (N-1)TS. 
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The  shape  of  the  power  spectral  density  of  the  desired  signal,  except  for  a 
gain  factor,  is  little  affected  by  the  nonlinearity.  This  is  illustrated  in  Figure  7.16  for 
the  case  where 

N=1 5,  L=10,  Og=5,  a=0.4975,  an=0.1,  Oc=5,  pc=0.9,  as=1.0,  fc/fs=0.5 
and  fs=1/Ts  is  the  pulse  repitition  frequency  of  the  transmitted  signal. 


Figure  7.16-  Power  spectral  densities  of  the  signal  at  the  input 
of  the  transversal  filter  when  RF  amplifier  is: 

(1)  linear ;  (2)  nonlinear. 

Because  the  desired  signal  power  spectral  density  takes  on  negative  values 

due  to  truncation  of  the  autocorrelation  sequence,  10Log10  Ssy(£2)  and 
-Nl 

10Log10  Ssj  (Q)  are  plotted  in  Fig.7.16.  Observe  that  the  output  power  spec¬ 
tral  density  is  identical  to  the  input  power  spectral  density  except  for  an  offset  of 
approximately  -7  dB.  This  result  is  readily  explained.  Eq.(7.57)  can  be  rewritten  as 
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With  reference  to  Eq.(6.12)  and  recalling  that  a=ag/o1(  the  above  expression  be¬ 
comes 

r  .2k +1  -  °2kl 


®  ■«  _i_  v  k  k 

(M  )  = - -  (M  )  1+  \  bk — -7- 

S  UV  1+a2  S  uv  k“1  k  22k 


1  l2k  as 


1  +a‘ 


(7.60) 

In  our  work  it  is  assumed  that  a2  «  a2.  Hence,  the  terms  involved  in  the  summa- 

s  c 

tion  on  k  are  negligible  and  we  obtain 


(Mf)  , 

5  uv  i  +  a2  5  uv 


(7.61) 


Consequently,  the  two  spectra  are  identical  except  for  a  multiplicative  constant. 
For  the  curves  shown  in  Fig.7.16,  a=0.4975.  The  -7  dB  offset  is  accounted  for  by 
the  fact  that 


lOLog 


=  -7.02dB  . 


a  =  0.4975 


(7.62) 


From  Eqs.(6.40)  and  (6.29),  the  uvth  element  of  mJJl  is  given  by 


,MNL.  ,mnl^  ,mNL. 
<Md  >„  =  (M  'uv-^s  >. 


<*  ..NL. 

o  )  nw 

1+a  normalized 


where 


(MNL) 


oo 

G2  Y  b, 

n2k  + 1 

°k 

1/  (2c2) 

normalized 

°  2-  Dk 

k  =  o 

22k  +  2 

1+a2 

2k- 

p  2k  -f  i 

it(v-u)TsJ| 

Ri  [(v- 

u)Ts]- 

<2os>  • 
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(7-63) 


(7.64) 


Comparison  of  Eqs.  (7.61)  and  (7.63)  reveals  that  o?/(  1+a2)  is  a  common 
multiplicative  factor  to  (m^l)  and  (m^l)  .  As  a  result,  the  signal-to-distur- 

s  uv  °  uv 

bance  ratio  given  by  Eq.  (6.39)  is  independent  of  this  factor. 

The  power  spectral  density  of  the  disturbance  at  the  nonlinearity  output  can 
be  expressed  as 

SdT(O)  =  i  (M"h  e-i(1-u>Q+  X  (MJJS  e-i,v-’>a  . 

ufN  d  u<  vT,  d 

(7-65) 

As  with  Eq.(7.58),  the  summations  over  u  and  v  incorporate  the  negative  and  non¬ 
negative  lags,  respectively. 

Both  Equations  (7.58)  and  (7.60)  contain  the  common  multiplicative  factor 
a2/(1  +a2).  In  order  to  compare  the  spectra  for  cases  II  and  III  with  those  for  case 
I,  it  is  convenient  to  divide  the  spectra  for  cases  II  and  III  by  the  factor  a2/(1+oc2). 
This  constrains  the  spectra  for  cases  II  and  III  at  the  input  to  the  transversal  filter 
to  have  identical  maxima  as  those  for  case  I.  Thereby,  any  differences  in  spectral 
shape  for  cases  II  and  III  from  those  of  case  I  are  readily  observed.  This  normal¬ 
ization  procedure  has  been  carried  out  in  all  of  the  remainig  spectral  plots.  With 
this  normalization,  the  desired  signal  spectrum  at  the  input  to  the  transversal  filter 
becomes  identical  for  cases  I,  II,  and  III. 

For  the  same  nonlinearity  as  in  Fig.  7.16,  Figs.  7.17  (a)  and  (b)  show  the  plot 
of  the  power  spectral  density  of  the  disturbance  at  the  input  to  the  transversal  filter 
before  and  after  normalization,  respectively.  When  normalized,  the  disturbance 
spectrum  for  case  II  has  the  same  maximum  as  the  disturbance  spectrum  for 
case  I. 

The  nonlinearity  tends  to  widen  the  power  spectral  density  of  the  distur¬ 
bance  (i.e,  clutter  plus  noise).  This  is  illustrated  in  Figure  7.17  (b).  Because  the 
contributions  of  the  nonlinear  terms  are  small  relative  to  the  linear  term,  the  output 
power  spectral  density  at  the  maximum,  in  the  vicinity  of  ft=0,  is  little  changed 
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Figure  7.17-  Power  spectral  densities  of  the  disturbance  at  the 


input  of  the  transversal  filter  when  RF  amplifier  is: 


' 1 )  linear,  (2)  nonlinear. 


a)  unnormalized  case,  ( b )  normalized  case. 


from  that  at  the  input.  However,  at  the  higher  frequencies  where  the  power  spec¬ 
tral  density  is  relatively  small,  the  output  power  spectral  density  is  in  excess  of  20 
dB  larger  than  that  at  the  input. 


7.5  Power  Spectral  Density  at  the  Transversal  Filter  Output 


Let  Sj  (Q)  denote  the  power  spectral  density  of  the  total  signal  at  the  out- 
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put  of  the  nonlinearity  (i.e.,  input  to  the  transversal  filter).  Because 


+  M 


NL 


(7.65) 


it  follows  that 


-  ML  -  IXL.  -  I  'I  L.  _ _ 

Sj  (£2)  =  Ssj  (£2)  +  Sdy  (£2)  .  (7.66) 

The  linear  transfer  function  of  the  transversal  filter  is  given  by  Eq.(7.13).  Let 
S0  (£2)  denote  the  power  spectral  density  of  the  transversal  filter  output.  Then 

2  -NL 

S0(£2)  —  JN-|-p(Q)|  •  Sy  (£2) 

2  -NL  -NL 

=  |HTF(£2)|2(SsT(£2)+SdT(£2))  .  (7.67) 

The  output  power  spectral  density  of  the  desired  signal  is 

o  -NL 

Sos(£2)  =  |Htf(Q)|2  SsT  (£2)  (7.68) 

whereas  the  power  spectral  density  of  the  disturbance  is 

2  -NL 

Sod(£2)  =  |HTF(£2)|2-SdT(£2)  -  (7.69) 

To  illustrate  the  effect  of  filtering  by  the  transversal  filter,  consider  the  same 
case  as  was  used  in  obtaining  the  plots  of  Figs.  7.16  and  7.17.  In  Fig.  7.18  plots 
are  shown  of  the  signal  and  disturbance  power  spectral  densities  at  the  filter  input 
plus  the  magnitude  of  the  filter  transfer  function  optimized  to  maximize  the  signal- 
to-disturbance  ratio  with  the  nonlinearity  taken  into  account .  The  filter  output  sig¬ 
nal  and  disturbance  power  spectral  ounsities  are  shown  in  Fig.  7.19.  Observe  that 
the  optimum  filter  attenuates  signals  at  low  frequencies  where  the  disturbance  is 
dominant  and  passes  signals  in  the  frequency  band  where  the  signal  spectrum  is 
maximum.  This  results  in  a  maximum  value  for  the  signal-to-disturbance  ratio. 
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Chapter  8 

COMPUTER  GENERATED  RESULTS 


In  order  to  assess  radar  performance  degradation  caused  by  the  receiver 
nonlinearity,  computer  generated  resuits  are  presented  and  discussed  in  this 
chapter.  Both  signal  and  receiver  parameters  are  varied  so  as  to  determine  how 
performance  is  affected.  In  addition  to  evaluating  signal-to-disturbance  ratios  at 
the  transversal  filter  output,  signal  spectra  and  transfer  functions  are  plotted 
where  ever  they  are  useful  in  interpreting  results.  In  terms  of  the  normalized  fre¬ 
quency  variable 


n  =  2rcfT  =  2n]~,  (8.1) 

fs 

the  spectra  and  transfer  functions  are  periodic  with  period  2jc.  Consequently,  they 
need  to  be  plotted  only  over  the  interval 


i  . 

Division  by  n  results  in  the  inequality 


-1  s 


f 


V2 


si 


Hence,  the  end  points  of  the  plots  correspond  to 

f  =  ±!?. 


The  variable 


is  used  in  all  of  the  plots. 


n  =  f 

*  V2 


(8-2) 


(8.3) 


(8-4) 


(8.5) 


99 


Three  different  cases  are  considered  depending  upon  whether  the  transver¬ 
sal  filter  weights  apply  to 

1)  an  ideal  linear  receiver  (case  I), 

2)  a  nonlinear  receiver  having  the  same  transversal  filter  as  was  used  with 
the  ideal  linear  receiver  (case  II), 

3)  a  nonlinear  receiver  with  the  transversal  filter  modified  to  account  for  the 
receiver  nonlinearity  (case  III). 

Unless  stated  otherwise,  a  situation  where  the  desired  signal  is  strong  rela¬ 
tive  to  the  thermal  noise  but  is  weak  relative  to  the  clutter  is  studied.  In  particular, 
the  standard  deviations  of  the  signal,  clutter,  and  noise  are  specified  to  be 


°s  =  '>ac  =  10’°n  =  01  •  (8.6) 

Consequently,  the  input  signal-to-noise  ratio  is 

(SNR)|n  =  20Log10^j  =  20dB  ,  (8.7) 

the  input  signal-to-clutter  ratio  is 

(SCR)in  =  20Log10(^j  =  -20dB  ,  (8.8) 

the  input  noise-to-clutter  ratio  is 

<NCR)in  =  20Log10^)  =  -40dB  ,  (8.9) 

and  the  input  signal-to-disturbance  ratio  is 

(SDR)jn  =  20Log10|~j  =  -20.044dB  (8.10) 


where 


(8.11) 


_2  _2  ,  _2  ,  2 
C1  =  as  +  °c+°n  " 


8.1  ComDuter  Generated  Results  for  Case  I:  The  Ideal  Linear  Receiver 


Case  I  refers  to  the  situation  where  the  radar  receiver  is  modeled  as  though 
it  is  linear  with  a  gain  of 


N*  o„ 


(8.12) 


Assuming  the  weights  of  the  transversal  filter  are  chosen  to  maximize  the  signal- 
to-disturbance  ratio  at  the  output  of  the  transversal  filter,  it  is  shown  in  Eq.(6.65) 
that  the  maximum  output  signal-to-disturbance  ratio  is  given  by 


G2a2 

(SDR)j  =  — ^  3T(M^)  s’ 


(8.13) 


where 


9  —  [a-,a2 . i 


(8-14) 


j2jcfn  (j-i)T 

=  e  D  s  ;  k=  1,2 . N 


(8-15) 


and  the  uv,h  element  of  the  matrix  MdL  is  given  in  Eq.(6.13)  as 


( (v-u)Tg)‘ 


<Md>, 


G2  2  a 
2  °c 


2  sin  (rcW  (v-  u)T$) 
7n  7iW(v-u)Te 


(8.16) 
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It  follows  that  each  element  in  is  proportional  to  2/G2.  Consequently, 

(SDR)j ,  where  the  subscript  I  refers  to  case  I,  is  independent  of  G  and,  therefore, 
is  also  independent  of  the  choices  for  L  and  og. 

Dependence  of  (SDR)j  on  the  remaining  parameters  is  difficult  to  ascertain 
by  inspection  because  of  the  complexity  of  the  (NxN)  matrix  (Mj1-)*1-  Therefore, 
computer  generated  results  are  used  in  the  following  subsections  to  study  the  ef¬ 
fects  of  these  parameters. 


8.1.1  Dependence  of  (SDR)r  on  Doppler  Frequency 

As  printed  out  in  chapter  7,  the  desired  signal  power  spectral  density  is  cen¬ 
tered  at  fo-  Because  one  period  of  the  spectra  spans  frequencies  in  the  interval 
0<  f<  fs,  frequencies  higher  than  fs  are  folded  over  into  this  interval.  To  avoid  am¬ 
biguities  in  the  Doppler  frequency,  it  is  assumed  that 


0  <^£1 


(8.17) 


A  plot  of  (SDR)!  versus  fp/fs  is  shown  in  Fig.8.1  for  three  different  values  of 


Pc 


(8.18) 


When  p^O,  the  clutter  power  spectral  density  is  white  (i.e.,  constant).  Since 
the  noise  power  spectral  density  is  also  white,  the  total  power  spectral  density  of 
the  disturbance  is  a  constant  with  respect  to  frequency.  Therefore,  (SDR)i  is  inde¬ 
pendent  of  the  center  frequency  of  the  desired  signal  spectrum  and  plots  as  a 
constant.  Even  so,  observe  that  the  transversal  filter  is  able  to  increase  the  output 
signal-to-disturbance  ratio  from  approximately  -20  dB  at  its  input  to  approximately 
-8  dB  at  its  output.  This  amounts  to  a  processing  gain  of  approximately  12  dB. 
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f(A 

Figure  8. 1  -  Dependence  of  (SDR),  on  the  Doppler  Frequency 

Wj>c=0 ,  (2) 

When  pc=0.9,  the  clutter  power  spectral  density  is  relatively  narrow.  Be¬ 
cause  of  periodicity  in  the  spectrum,  it  peaks  both  at  fo/fs=0  and  fo/fs=1.  The 
smallest  values  of  the  power  spectral  density  occur  in  the  neighborhood  of 
fs=0.5.  Consequently,  the  plot  of  (SDR)!  is  maximum  at  this  value  of  Doppler  fre¬ 
quency.  When  fQ=0.5fs,  a  processing  gain  in  excess  of  50  dB  is  achieved.  (Note 
that  the  input  signal-to-disturbance  ratio  has  been  increased  from  approximately 
-20  dB  to  approximately  31  dB). 

pc=0.5  corresponds  to  an  intermediate  situation.  The  smallest  values  of  the 
clutter  power  spectral  density  are  still  in  the  vicinity  of  fo/fs=0.5  so  that  the  plot  of 
(SDR)!  peaks  at  this  value.  However,  because  the  clutter  power  spectral  density  is 
much  wider  than  it  was  for  pc=0.9,  the  clutter  spectrum  is  significantly  larger  at 
fo=0.5fs  and  the  processing  gain  is  only  approximately  20  dB. 

8.1.2  Dependence  of  (SDR)j  on  The  Clutter  Correlation  Parameter 

As  pointed  out  in  chapter  7,  the  clutter  correlation  sequence  can  be  ex- 
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pressed  as 


Rc[kTs]  =  2o|-pf  (8.19) 

where  pc  is  given  by  Eq.(8.18).  The  clutter  samples  are  said  to  be  highly  correlat¬ 
ed  when  pc=1  and  uncorrelated  when  pc=0.  A  plot  of  the  normalized  clutter  corre¬ 
lation  sequence,  (Rc  [kTs))  /  (2o^) ,  is  shown  in  Fig.8.2.  The  correlation 
sequence  is  observed  to  become  wider  as  the  value  of  pc  is  increased.  As  a  re¬ 
sult,  the  clutter  power  spectral  density  becomes  narrower  for  larger  values  of  pc.  A 
plot  of  (SDR)!  versus  pc  is  shown  in  Fig. 8. 3  for  three  different  values  of  fo/fs. 


Figure  8.2  -  Normalized  Clutter  Correlation  Sequence 


As  shown  in  Fig.8.1,  the  most  favorable  positioning  of  the  desired  signal 
spectrum  occurs  when  fo=0.5fs.  Because  the  clutter  spectrum  becomes  narrower 
as  pc  becomes  larger,  the  signal-to-disturbance  ratio  increases  with  increasing 
values  of  pc. 
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Since  the  clutter  spectrum  peaks  at  f=0,  the  worst  case  positioning  of  the  de¬ 
sired  signal  spectrum  occurs  when  fD=0.  The  area  of  the  clutter  power  spectral 
density  remains  constant  at  2oc2  independent  of  the  value  for  pc.  Therefore,  as 
the  clutter  spectrum  becomes  narrower,  its  maximum  value  at  f=0  becomes  larger. 
As  a  result,  the  output  signal-to-disturbance  ratio  decreases  as  pc  increases. 

An  intermediate  situation  results  when  fD=0.25fs.  Because  the  peak  of  the 
desired  signal  spectrum  is  offset  from  that  of  the  clutter  spectrum,  the  output  sig¬ 
nal-to-disturbance  ratio  improves  as  the  clutter  spectrum  narrows.  Nevertheless, 
performance  is  not  as  good  as  when  fD=0.5fs. 

For  pc=0,  the  clutter  power  spectral  density  is  constant.  Consequently,  the 
location  of  the  desired  signal  power  spectral  density  is  immaterial  and  all  three 
curves  in  Fig.8.3  have  identical  values  at  pc=0. 
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From  Eq.(7.29),  the  correlation  sequence  of  the  thermal  noise  is 


sin  (rcWkT  ) 
Rn[KTs]=2^— 


As  explained  in  Chapter  7,  this  can  be  approximated  by 


(8.20) 


{2  c2  ;  k=0 

"  (8.21) 

0  ;  kiO  - 

As  a  result,  the  thermal  noise  power  spectral  density  is  approximately  constant. 
Throughout  this  work,  it  is  assumed  that  the  total  clutter  power  is  much  larger  than 
the  total  thermal  noise  power.  Consequently,  as  cn  is  varied  over  the  interval 

0.1$onS1  ,  (8.22) 

cc  is  fixed  at  cc=10. 

Because  the  clutter  power  spectral  density  is  Gaussian  shaped  while  the 
noise  power  spectral  density  is  constant,  it  is  possible  for  the  thermal  noise  to  be 
significant  when  the  desired  signal  spectrum  is  positioned  at  the  minimum  of  the 
clutter  power  spectral  density.  This  is  illustrated  in  Fig. 8.4  where  pc=0.9  and 
fo=0.5fs.  The  power  spectral  densities  at  the  input  to  the  transversal  filter  are- 
shown  in  parts  (a)  and  (b)  of  the  figure  while  the  output  power  spectral  densities 


Figure  8.4  -  a)  Power  spectral  density  of  desired  signal  at  input  to  the 
transversal  filter  ( fn=0.5L) .  b)  power  spectral  density  of  clutter  at  input 
to  the  transversaffilter  fb^=0.9).  c)  power  spectral  densities  at  output 
of  transversal  filter  foron=0. 1:  (1 )  clutter  plus  noise.  (2)  desired  signal, 
(d)  power  spectral  densities  at  output  of  transversal  filter  foro„=1: 

(1)  clutter  plus  noise.  (2)  desired  signal. 
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are  shown  in  parts  (c)  and  (d)  for  the  noise-to-clutter  ratio  equal  to  -40dB  and  -20 
dB,  respectively.  Note  that  the  transversal  filter  is  successful  in  removing  most  of 
the  low-frequency  clutter  power.  However,  the  increase  in  thermal  noise  power  of 
20dB  does  result  in  an  increase  of  the  noise  floor  by  approximately  20dB  at  the 
transversal  filter  output.  For  an=0.1,  (SDRJpSO.BdB  while  (SDR)!=11.4dB  for 
crn=1 .  Hence,  the  20dB  increase  in  thermal  noise  power  results  in  a  degradation 
of  19.2dB  in  the  output  signal-to-disturbance  ratio. 

On  the  other  hand,  when  the  desired  signal  spectrum  is  positioned  at  the 
peak  of  the  clutter  power  spectral  density,  the  20dB  increase  in  thermal  noise 
power  has  negligible  effect.  This  is  illustrated  in  Fig.8.5  where,  as  in  Fig. 8.4, 
pc=0.9  but  now  fD=0.  Observe  that  the  output  power  spectral  densities  in  part  (d) 
are  little  changed  from  those  in  part  (c).  For  an=0.1,  (SDR)p-14.2dB  while 
(SDR)j=-14.5dB  for  on=1 .  Thus,  the  20dB  increase  in  thermal  noise  power  results 
in  only  a  0.3dB  decrease  of  the  output  signal-to-disturbance  ratio. 

Analogous  to  the  results  of  Fig.8.5,  when  fD=0,  the  20dB  increase  in  thermal 
noise  power  has  negligible  effect  when  the  clutter  power  spectral  density  is  white. 
This  is  illustrated  in  Fig. 8.6  where  pc=0  and  fo=0.25fs.  As  in  Fig.8.5,  the  output 
power  spectral  densities  in  part  (d)  are  little  changed  from  those  in  part  (c).  For 
an=0.1,  (SDR)!=-8.2dB  while  (SDR)p-8.3dB  for  an=1.  Here,  the  20dB  increase  in 
thermal  noise  power  results  in  only  a  0.1  dB  decrease  of  the  output  signal-to-dis- 
turbance  ratio. 
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Figure  8.5  a)  Power  spectral  density  of  desired  signal  at  input  to  the 
transversal  filter  (fn=0.5fJ.  b)  power  spectral  density  of  clutter  at  input 
tg  the  transversaffilter  (o~=0.9).  c)  power  spectral  densities  at  output 
of  transversal  filter  fore 2=0. 1:  (1)  clutter  plus  noise.  (2)  desired  signal . 
(d)  power  spectral  densities  at  output  of  transversal  filter  for  o2=1: 

(1)  clutter  plus  noise.  (2)  desired  signal. 
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Plots  of  (SDR){  versus  the  input  noise-to-clutter  ratio,  as  defined  in  Eq.(8.9), 
are  shown  in  Fig. 8. 7.  In  part  (a),  fD=0.5fs.  For  highly  correlated  samples  (pc=0.9), 
where  the  clutter  power  spectral  density  is  relatively  narrow,  the  output  signal-to- 
disturbance  ratio  falls  off  linearly  in  dB  as  the  thermal  noise  power  is  increased. 
Specifically,  a  20  dB  increase  in  (NCR)in  results  in  approximately  a  20  dB  de¬ 
crease  of  (SDR)j.  For  moderately  correlated  (pc=0.5)  and  uncorrelated  (pc=0) 
clutter,  where  the  clutter  power  spectral  density  is  wider,  the  signal-to-disturbance 
ratio  remains  constant  as  an  is  increased.  However,  (SDR)j  is  approximately  8  dB 
larger  at  pc=0.5  than  it  is  at  pc=0.  In  part  (b)  of  the  figure,  fD=0.25fs.  For 
pc=0.9,the  decrease  in  (SDR)j  is  approximately  3.5  dB, 


CD 

V 


(d) 
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Figure  8.6  a)  Power  spectral  density  of  desired  signal  at  input  to  the 
transversal  filter  (fn=0.25L).  b)  power  spectral  density  of  clutter  at  in¬ 
put  to  the  transversal  filterlo^=0).  c)  power  spectral  densities  at  output 
of  transversal  filter  for  <3  n=0.  T:  (1)  clutter  plus  noise.  (2)  desired  signal, 
(d)  power  spectral  densities  at  output  of  transversal  filter  for  1 : 

( 1 )  clutter  plus  noise.  (2)  desired  signal. 

considerably  less  than  it  was  for  fD=0.5  fs.  Also,  the  constant  output  signal-to-dis- 
turbance  ratio  at  pc=0.5  is  only  approximately  1  dB  larger  than  the  value  for  pc=0. 
In  part  (c)  of  the  figure,  fo=0.  Thus,  the  clutter  and  desired  signal  power  spectral 
densities  both  peak  at  f=0.  For  pc=0.9,  (SDR)j  is  relatively  insensitive  to  the  in¬ 
crease  in  thermal  noise  power  (note  the  expanded  vertical  scale)  because  the  dis¬ 
turbance  is  dominated  by  the  clutter  in  the  vicinity  of  f=0.  As  in  parts  (a)  and  (b), 
the  curves  are  constant  for  pc=0.5  and  pc=0.  However,  in  contrast  to  parts  (a)  and 
(b),  the  output  signal-to-disturbance  ratio  decreases  with  increasing  values  of  pc. 
As  explained  previously,  this  is  due  to  the  fact  that  the  peak  value  of  the  clutter 
power  spectral  density  increases  at  f=0  with  increasing  values  of  pc  because  the 
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(SDR),-dB-  (SDR),-dB- 


area  under  the  density  remains  constant  as  the  spectrum  becomes  narrower. 
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(C) 


Figure  8. 7  -  Dependence  of  {SDR),  on  Thermal  Noise  Power  Level: 

(1)  o ,=0.9,  (2)  o >0.5,  (3)  o,=0 
(a)  fn/L=0.5,  (b)  fn/M.25,  ( c)ln/t c^g 

8.1.4  Dependence  of  (SDR)t  on  the  Clutter  Power  Level 

In  this  section  oc  is  varied  over  the  interval 

10  <  ac  <  100  (8.23) 

while  os  and  on  are  held  fixed  at  as=1  and  cn=0.1 .  Therefore,  the  input  signal-to- 
clutter  ratio  and  input  noise-to-clutter  ratio,  as  defined  by  Eqs.  (8.8)  and  (8.9),  vary 
according  to 


-40dB  <  (SCR)jn  <-20dB  (8.24) 

and 

-60dB^  (NCR)jn^-40dB  .  (8.25) 

Dependence  of  (SDR)j  on  the  clutter  power  level  is  illustrated  in  Fig.8.8.  In 
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part  (a)  of  the  figure,  fD=0.5fs.  For  pc=0.9,  the  clutter  power  spectral  density  is  suf¬ 
ficiently  narrow,  such  that  the  thermal  noise  dominates  the  disturbance  at  f=0.5fs. 
Consequently,  the  output  signal-to-disturbance  ratio  is  relatively  insensitive  to 
variations  in  the  clutter  power.  However,  for  the  remaining  cases  considered  in 
parts  (a),  (b),  and  (c),  the  clutter  dominates  the  disturbance  in  the  frequency  inter- 
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(C) 


Figure  8.8  -  Dependence  of  (SDR),  on  the  Clutter  Power  Level: 
jjggsL-g n-0.1):  (D  o^O.Q; (2)  o^0.5,  (3)  o^O 
-  (31  frr^O.5.  (b)lr/fz=0.25. 70  Ir/f^O.  ~ 

val  where  the  major  part  of  the  desired  signal  power  is  located.  As  a  result,  the 
curves  plot  as  straight  lines  with  slopes  such  that  a  20  dB  increase  in  the  input  sig- 
nal-to-clutter  ratio  is  accompanied  by  approximately  a  20  dB  increase  in  the  out¬ 
put  signal-to-disturbance  ratio.  This  is  reasonable  because  (SDR)!  is  given  by 
Eq.(8.13)  where  the  uvth  element  of  the  matrix  Mj-  is  given  by  Eq.(8.16).  For 
those  frequencies  where  the  clutter  power  is  much  larger  then  the  thermal  noise 
power,  the  matrix  (MdL)_1  is  inversely  proportional  to  ac2. 

8.1.5  Dependence  of  (SDR)t  on  the  Number  of  Samples.  N. 

It  is  expected  that  the  output  signal-to-disturbance  ratio  should  increase  with 
an  increase  in  the  number  of  samples  processed.  This  is  shown  to  be  the  case  in 
Fig.  8.9  for  fo/fs=0.5,  0.25  ,0  and  pc=0.9,  0.5,  0. 

The  improvement  is  due  to  two  factors:  1 )  the  power  in  the  desired  signal  be¬ 
comes  more  concentrated  in  the  frequency  interval  centered  at  fp,  and  2)  the 
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number  of  degrees  of  freedom  of  the  transversal  filter  is  equal  to  N.  The  input 
power  spectral  densities  for  N=15  and  N=35  are  shown  in  Fig.  8.10.  Note  that  the 
clutter  plus  noise  power  spectral  density  remains  unchanged  whereas  the  power 
spectral  density  of  the  desired  signal  has  a  narrower  central  lobe  with  its  peak 
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value  having  increased  by  approximately  3  dB.  The  corresponding  output  power 
spectral  densities  are  shown  in  Fig.  8.11.  The  increased  degrees  of  freedom  of 
the  transversal  filter  result  in  spectra  having  finer  grain  where  the  clutter  plus 
noise  density  has  been  reduced  approximately  1 0dB  in  the  vicinity  of  f=0  and  3  dB 
in  the  vicinity  of  f=0.5fs. 

8.2  Computer  Generated  Results  for  Case  II;  The  Nonlinear  Receiver 
with  the  Same  TVansversal  Filter  as  in  Case  I 

Case  II  refers  to  the  situation  where  the  radar  receiver  is  nonlinear.  Neverth- 
less,  the  transversal  filter  designed  for  the  ideal  linear  receiver  is  used  to  improve 
the  output  signal-to-disturbance  ratio.  This  case  can  arise  when  either  the  design 
engineer  is  unaware  of  the  nonlinearity  or  does  not  choose  to  take  it  into  account. 

From  Eq.(6.67)  the  output  signal-to-disturbance  ratio  is  given  by 
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(SDR) 


WHMg  LW 

11  wHm^lw 


(8.26) 


where  the  transversal  filter  weight  vector  is  determined  according  to  Eq.(6.63)  as 

-i 

(8.27) 


W  =  k-  (Mjp  s*  . 


In  Eq.(8.27)  k  is  an  arbitrary  constant  and  3  is  defined  by  Eqs.(8.14)  and  (8.15). 
For  simplicity,  k  is  set  equal  to  unity.  M^1-  is  the  (NxN)  matrix  whose  uvth  element 
is  given  by  Eq.(8.16).  The  uvth  element  of  the  (NxN)  matrix  MSNL  is  given  by 
Eq.(6.29)  as 
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00 


p2k  + 1 

S _ 

22k+i 
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Finally,  the  uvth  element  of  the  (NxN)  matrix 
(6.40)  as 


-ei2n'D""U,Ts  .  (8.28) 

NL  is  given  by  Eqs.  (6.34)  and 


<Mj/L)  =  (MNL)  -(M^L)  (8.29) 

u  UV  uv  O 
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and,  from  Eq.(6.25), 


In  Chapter  6  it  is  shown  that  suitable  approximations  to  the  infinite  summations  in 
Eqs.  (8.28)  and  (8.30)  can  be  achieved  by  summing  over  the  first  12  terms. 

By  substitution  of  Eq.(8.27)  into  Eq.(8.26),  the  expression  for  the  output  sig- 
nal-to-disturbance  ratio  becomes 


(SDR)  I]t  = 


“1  —1 

aT(M|j)  s 

__  __  . 

ST<mL)  M^L(mL)  a 


(8.32) 


Because  all  matrices  are  proportional  to  G2,  (SDR)U  is  independent  of  G.  From 
Eq.(8.12),  it  follows  that  the  output  signal-to-disturbance  ratio  is  independent  of 
the  saturation  level,  L,  of  the  nonlinearity. 


8.2.1  Dependence  of  (SDR)jt  on  the  Clutter  Correlation  Parameter 

As  pointed  out  in  Chapter  7,  when  the  clutter  is  much  larger  than  the  signal, 
the  nonlinearity  tends  to  broaden  the  clutter  power  spectral  density  without  signifi¬ 
cantly  changing  the  shape  of  the  desired  signal.  Tnis  is  illustrated  in  Fig.8.12  for 
the  situation  where 

os=1 ,  ac=10,  on=0.1,  cg=5  ,  asOg/a!  =0.4975. 

Since  the  standard  deviation  of  the  total  input  signal  is  approximately  twice  the  lin¬ 
ear  range  of  the  nonlinearity,  the  input-output  characteristic  of  the  nonlinearity  is 
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Figure  8. 12  -  Power  Spectral  Densities  at  Nonlinear  Output 
(o«=1.  02=5.  a=0.4975):(1)  easel  .(2)  case II 
(a)  Desired  signal  (fn/L=0.5).  (b)  Desired  signal  (fp/fz=0.25). 

(c)  clutter  plus  noise  (d)  clutter  plus  noise  (oz=0). 

being  driven  into  its  nonlinear  region.  However,  note  that  the  desired  signal  is 
small  enough  such  that,  by  itself,  it  would  operate  entirely  in  the  linear  region  of 
the  nonlinearity.  In  part  (a)  where  fo=0.5fs  and  in  part  (b)  where  fD=0.25fs,  the 
shape  of  the  power  spectral  density  of  the  desired  signal  at  the  nonlinearity  output 
is  seen  to  be  unchanged.  This  is  due  to  the  relatively  small  desired  signal  level.  In 
part  (c),  where  pc=0.9,  the  power  spectral  density  of  the  clutter  plus  noise  at  the 
nonlinearity  output  is  observed  to  be  noticeably  wider  due  to  the  nonlinearity.  Be¬ 
cause  the  nonlinear  terms  in  Eq.(8.30)  are  much  smaller  than  the  linear  term,  the 
most  noticeable  increase  in  the  spectral  level  occurs  at  high  frequencies  where 
the  disturbance  power  spectral  density  is  small.  At  f=fs,  the  increase  is  approxi¬ 
mately  22  dB  whereas  at  f=0,  it  is  on  the  order  of  a  few  tenths  of  a  dB .  In  part  (d), 
where  pc=0  and  the  disturbance  power  spectral  density  is  a  constant  both  at  the 
input  and  output  of  the  nonlinearity,  the  increase  in  spectral  level  is  seen  to  be 
around  0.4  dB. 

Dependence  of  the  output  signal-to-disturbance  ratio  on  the  clutter  correla¬ 
tion  parameter,  pc,  is  shown  for  both  case  I  and  case  II  in  Fig.8.13.  In  part  (a)  of 
the  figure,  fo=0.5fs.  Here  the  desired  signal  spectrum  is  located  at  the  minimum 
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value  of  the  clutter  plus  noise  power  spectral  density.  Consistent  with  the  in¬ 
creased  levels  of  the  disturbance  seen  in  Fig.8.12,  the  nonlinearity  causes  a  22 
dB  loss  of  the  output  signal-to-disturbance  ratio  for  pc=0.9  and  an  approximate  0.4 
dB  loss  for  pc=0.  This  is  to  be  expected  since  the  transfer  function  of  the  transver¬ 
sal  filter  is  the  same  for  both  cases  I  and  II.  In  Fig.8.13  (b),  where  fD=0.25fs,  the 
losses  are  approximately  10  dB  for  pc=0.9  and  0.4  dB  for  pc=0.  The  10  dB  Joss  is 
consistent  with  the  increase  in  disturbance  level  shown  in  Fig.8.12  (c)  at  12/7t=0.5 
(i.e.;  f=0.25fs).  In  Fig.8.13  (c),  fD=0.  The  decrease  in  (SDR)u  for  increasing  values 
of  pc  follows  the  same  trend  observed  in  the  linear  case.  However,  the  loss  of  ap¬ 
proximately  2.2  dB  for  pc=0.9  is  larger  than  the  0.4  dB  loss  predicted  at  f=0  in 
Fig.8.12  (c).  Note  that 

(SDR)j>  (SDR)n  (8-33) 

for  all  values  of  pc  in  Figs.  8.13  (a),  (b),  and  (c).  This  is  reasonable  since  the 
shape  of  the  power  spectral  density  of  the  desired  signal  has  been  unaffected  by 
the  nonlinearity  while  the  width  of  the  disturbance  power  spectral  density  has 
been  increased. 


(a) 


Figure  8.13-  ( Caption  appears  on  next  page). 
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Recall  that  cg  is  a  measure  of  the  linear  range  of  the  nonlinearity  whereas  a ( 
represents  the  r.m.s.  value  of  the  total  input  signa*.  Consequently,  the  larger  is  the 
value  of  a,  the  more  linear  is  the  behavior  of  the  nonlinearity.  As  a  increases,  it  is 
expected  that  the  output  signal-to-disturbance  ratio  for  Case  II  will  approach  that 
of  Case  I.  This  is  verified  in  Fig.8.14  where  a  is  varied  by  changing  og  from  5  to  50 
while  all  signal  levels  are  held  fixed.  In  Fig.8.14  (a),  where  fo=0.5fs  and  pc=0.9, 
the  system  is  seen  to  behave  linearly  for  a>5.  For  all  of  the  other  cases  consid¬ 
ered  in  the  figure,  linearity  is  approached  for  a>3.  Also,  the  inequality  in  Eq.(8.33) 
is  satisfied  for  each  case. 


(a) 


(b) 


a 

Figure  8.14-  ( Caption  appears  on  page  128). 
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(SDR)i  -dB-  (SDR);  -dB-  (SDR)j  -dB- 
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As  in  section  8.1 .3,  it  is  assumed  that  the  total  clutter  power  is  always  much 
larger  than  the  total  thermal  noise  power.  Consequently,  ac  is  held  fixed  at  ac=1 0 
while  an  is  varied  over  the  interval 


0.1  <o  £1  . 
n 


(8.34) 


Plots  of  (SDR)n  versus  the  input  noise-to-clutter  ratio,  as  defined  in  Eq.(8.9), 
are  shown  in  Fig. 8.1 5.  To  provide  a  basis  for  comparison,  curves  of  (SDR)T  for  the 
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linear  case  are  also  shown.  Because  of  the  broadening  of  the  clutter  power  spec¬ 
tral  density  by  the  nonlinearity,  the  output  signal-to-disturbance  ratio  is  less  sensi¬ 
tive  to  an  increase  in  the  thermal  noise  power  than  was  the  case  for  the  linear 
receiver.  This  is  seen  in  Figs.8,15  (a)  and  (c).  In  some  cases,  the  output  signal-to- 
disturbance  ratio  actually  increases  with  an  increase  in  thermal  noise  power,  as 
seen  in  Figs.8.1 5  (c)  and  (e).  This  latter  behavior  can  be  understood  by  examining 
Fig.8.16  which  shows  some  of  the  input  spectra  and  the  corresponding  transfer 
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functions  related  to  the  plot  of  Fig.8.15  (e).  The  power  spectral  density  of  the  de¬ 
sired  signal  is  shown  in  Fig.8.16  (a) .  The  dutter  plus  noise  power  spectral  densi¬ 
ties  for  on=0.1  and  on=1  are  shown  in  Fig.8.16  (b)  and  (d),  respectively.  Finally, 
the  squared  magnitudes  of  the  transversal  filter  transfer  functions  for  cn=0.1  and 
an*1  are  shown  in  Figs.8.16  (c)  and  (e),  respectively.  When  on=0.1,  the  transver¬ 
sal  filter  tries  to  maximize  the  output  signal-to-disturbance  ratio  by  passing  fre¬ 
quencies  limited  to  the  range 


0.53fs^|f|<fs  (8.35) 

even  though  the  main  lobe  of  the  desired  signal  is  centered  at  fo=0.  When  an=1, 
the  tails  of  the  clutter  plus  noise  power  spectral  density  have  been  increased  by 
approximately  6  dB.  The  transversal  filter  is  now  able  to  maximize  the  output  sig¬ 
nal-to-disturbance  ratio  by  passing  frequencies  over  the  entire  frequency  range 
given  by 


0<|f|£fs.  (8.36) 

This  results  in  a  larger  value  for  (SDR)n  when  an=1  than  when  an=0.1 . 


131 


Figure  8.16-  Input  spectra  and  transfer  functions  related  to 
Fia.8. 15(e):  (a)  power  spectral  density  of  the  desired  signal .  lb)  clutter 
>. lus  noise  power  spectral  density  when  gn=0. 1.  (c)  squared  magnitude 
of  transfer  function  when  on=0. 1.  (d)  clutter  plus  noise  power  spectral 
density  when  o„=1.  (e)  squared  magnitude  of  transfer  function  when 

On  =7 

(1 )  easel.  (2)  case ll. 


As  in  section  8.1.4,  oc  is  varied  over  the  interval 


10  <  oc  <  100 

while  os  and  on  are  held  fixed  at  as=1  and  an=0.1 .  Dependence  of  (SDR)n  on  the 
clutter  power  level  is  illustrated  in  Fig.8.17.  To  provide  a  basis  for  comparison, 
curves  of  (SDR)j  for  the  linear  case  are  also  shown.  As  expected,  when  the  clutter 
dominates  the  disturbance  in  the  frequency  region  about  fD,  the  curves  decrease 
significantly  as  the  total  clutter  power  increases  by  20  dB.  As  oc  varies  from  10  to 
100,  a  varies  from  0.4975  to  0.05.  Thus,  the  nonlinearity  is 
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(SDR)i  -dB-  (SDR);  -dB-  (SDR);  -dB- 


(SDR)i-dB-  (SDR)i-dB-  (SDRfc-dB- 


driven  harder  as  the  clutter  power  is  increased.  This,  in  turn,  results  in  increased 
spreading  of  the  clutter  power  spectral  density.  Consequently,  compared  to  the  lin¬ 
ear  receiver,  the  degradation  in  output  signal-to-disturbance  ratio  is  greater  for  the 
nonlinear  receiver  (see  Figs.8.17  (a),  (c),  and  (e)  where  pc=0.9).  This  effect  is  not 
seen  in  Figs.8.17(b),  (d),  and  (f)  where  pc=0  and  the  disturbance  power  spectral 
density  is  white. 

8.2.5  Dependence  of  (SDR)tt  on  the  Doppler  Frequency 

A  plot  of  (SDR)u  versus  fo/fs  is  shown  in  Fig.8.18  for  pc=0.9  and  pc=0.  To 
provide  a  basis  for  comparison,  curves  of  (SDR)j  for  the  linear  case  are  also 
shown.  When  pc=0.9  both  (SDR)j  and  (SDR)n  peak  at  fo/fs=0.5.  This  corresponds 
to  the  situation  where  the  desired  signal  spectrum  is  positioned  at  the  minimum  of 
the  clutter  power  spectral  density.  When  pc=0,  the  clutter  plus  noise  power  spec¬ 
tral  density  is  white.  Hence,  the  output  signal-to-disturbance  ratios  are  indepen¬ 
dent  of  fo/fs-  in  all  cases  note  that  the  nonlinearity  results  in 

(SDR) j>  (SDR)jj  .  (8.37) 
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8.2.6  Dependence  Of  (SDR)™  on  the  Number  of  Samples,  N 

The  dependence  of  (SDR)u  on  the  number  of  samples,  N,  is  shown  in 
Fig.8.19  for  ^=0.5,  0.25,  0  and  pc=0.9,0.  Plots  for  (SDR)j  are  also  included  as  a 
basis  for  comparison.  As  expected,  the  output  signal-to-disturbance  ratios  in¬ 
crease  with  an  increase  in  the  number  of  samples  processed.  The  amount  of  im¬ 
provement  for  the  nonlinear  receiver  approximately  equais  that  of  the  linear  recei- 
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(SDR),  -dB- 
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ver.  As  explained  in  section  8.1 .5,  the  improvement  is  due  to  1 )  the  power  spectral 
density  of  the  desired  signal  becoming  more  concentrated  about  fD  and  2)  the  de¬ 
grees  of  freedom  of  the  transversal  filter  increasing  with  N. 

8.3  Computer  Generated  Results  for  Case  III:  Nonlinear  Receiver  with 


the  Transversal  Filter  optimized  to  Account  for  the  Receiver  Non 


linearity 


Case  in  refers  to  the  situation  where  the  radar  receiver  is  nonlinear  and  the 
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transversal  filter  is  optimized  to  account  for  the  receiver  nonlinearity  so  as  maxi¬ 
mize  the  output  signal-to-disturbance  ratio. 


From  Eq.(6.68)  the  output  signal-to-disturbance  ratio  is  given  by 


(SDR) 


wVLw 

O 

Ilf  "  -H..NL- 
W 


(8.38) 


where  the  transversal  filter  weight  vector  is  determined  according  to  Eq.(6.74)  as 


NL  "1  * 


W  =  k  •  (M^  )  3 


(8.39) 


In  Eq.(8.39)  k  is  an  arbitrary  constant  and  3  is  defined  by  Eqs.(8.14)  and  (8.15). 
For  simplicity,  k  is  set  equal  to  unity.  MSNL  is  the  (NxN)  matrix  whose  uvth  element 
is  given  by  Eq.(8.28).  The  uvth  element  of  the  (NxN)  matrix  MdNL  is  given  by 
Eq.(8.30).  These  elements  are  expressed  as  infinite  series.  Once  again,  suitable 
approximations  are  obtained  by  summing  over  the  first  12  terms  in  the  series. 

From  Eqs.(6.70)  and  (6.75),  the  expression  for  the  output  signal-to-distur¬ 
bance  ratio  can  be  simplified  to 

(SDR)m  =  A-31  (Mjj)  3  (8.40) 


where 


A 


r2k  + 1 
°k 


22k  +  i 


1/(2ap 

1  +a2 


n2k  +  i 


<2a‘) 


2k  +  i 


(8.41) 


Because  (M^1-)’1  is  inversely  proportional  to  G2  whereas  A  is  proportional  to  G2, 
(SDR)m  is  independent  of  G.  From  Eq.(8.12),  it  follows  that  the  output  signal-to- 
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disturbance  ratio  is  independent  of  L,  the  saturation  level  of  the  nonlinearity. 


8.3.1  Dependence  of  (SDR)ttt  on  the  Clutter  Correlation  Parameter 

Dependence  of  the  output  signal-to-disturbance  ratio  on  the  clutter  correla¬ 
tion  parameter,  pc>  is  shown  in  Fig. 8.20  for  cases  I,  II,  III.  In  part  (a),  where 
fo=0.5fs,  case  ID  shows  negligible  improvement  over  case  II  for  0<pc<0.6.  There¬ 
after,  the  improvement  increases  with  increasing  values  of  pc.  However,  at  Pc=0.9, 
the  improvement  is  only  around  1  dB.  In  parts  (b)  and  (c),  where  fQ/fs=»0.25  and  0, 
respectively,  case  III  provides  a  noticeable  improvement  over  case  II  for  pc>0.53 
which  increases  to  approximately  2  dB  at  pc=0.9.  It  is  concluded  that  the  more 
correlated  is  the  clutter,  the  better  is  case  III  able  to  outperform  case  II.  To  put  it 
another  way,  the  narrower  is  the  shape  of  the  clutter  power  spectral  density,  the 
more  distorted  is  the  spectrum  by  the  nonlinearity  and  the  more  that  can  be  done 
by  shaping  the  transfer  function  of  the  transversal  filter  to  maximize  the  output  sig- 
nal-to  disturbance  ratio. 


Fiaure  8.20  -  (Caption  appears  on  next  page 
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Figure  8.20  -  Dependence  of  Output  Sianal-to-Disturbance  Ratios  on 


the  Clutter  Correlation  Parameter  (oc-1 ,  a„=5.  a =0.4975 


B 


8.3.2  Dependence  of  (SDR)m  on  01=0^0^ 

As  explained  in  section  8.2.2,  the  larger  is  the  value  of  a,  the  more  linear  is 
the  behavior  of  the  nonlinearity.  Dependence  of  the  output  signal-to-disturbance 
ratio  on  a  is  shown  in  Fig.8.21  for  cases  I,  II,  and  III.  The  simulation  was  carried 
out  by  allowing  cg  to  vary  in  the  interval  5<og<50  while  all  signal  levels  were  held 
constant.  As  suggested  by  the  results  presented  in  Fig.8.20,  there  is  little  differ¬ 
ence  between  cases  II  and  III  when  pc=0.  For  pc=0.9,  case  III  is  better  able  to 
outperform  case  n  as  the  nonlinearity  is  driven  harder.  Consistent  with  Fig.8.20, 
the  improvement  for  pc=0.9  and  a=0.5  is  approximately  1  dB  for  fD*0.5fs  and  ap¬ 
proximately  2  dB  for  fo=0.25fs  and  0. 
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Figure  8.21  -  Dependence  of  (SDR),,,  on  a=a!Jc>,: 

(1)  (SDR),  -dB-.  (2)  (SDR),,  -dB-,  (3T(SDRh„  rdB-_ 
a)  Uf^O.5.  o^Q.9,  b)  WU0.5.  o ,=0  c)  M*=025,  p^0.9. 

(d)  lr/L=0.25.  pr=0.  H  Mz=0,  p >0.9, 7)  fn/f^O,  Qc=0. 

8.3.3  Dependence  of  (SDR)m  on  the  Thermal  Noise  Power  Level 

As  in  sections  8.1 .3  and  8.2.3,  it  is  assumed  that  the  total  clutter  power  is  al¬ 
ways  much  larger  than  the  total  thermal  noise  power.  Consequently,  oc  is  held 
fixed  at  oc=10  while  cn  is  varied  over  the  interval 

0.1<on<1.  (8.42) 

Plots  of  (SDR)j  versus  the  input  noise-to-clutter  ratio,  as  defined  in  Eq.(8.9), 
are  plotted  in  Fig.8.22  for  i=I,  II,  III.  Once  again,  there  is  little  difference  between 
cases  II  and  m  for  pc=0.  However,  for  pc=0.9,  case  III  performs  better  than  case 
II  as  the  thermal  noise  power  is  reduced.  Reduction  of  the  thermal  noise  power 
results  in  the  clutter  plus  noise  power  spectral  density  becoming  less  flat  in  the 
tails  of  the  spectrum.  The  plots  in  Figs.8.22  (c)  and  (e)  suggest  that  the  transver¬ 
sal  filter  transfer  function,  optimized  for  case  ni,  is  able  to  take  advantage  of  this 
change  in  spectral  shape. 
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(SDR)i  -dB-  (SDR)i  -dB-  (SDR)t  -dB- 
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(SDR)j  -dB- 
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8.3.4  Dependence  of  (SDR)m  on  the  Clutter  Power  Level 

As  in  sections  8.1.4  and  8.2.4,  as  and  an  are  held  constant  at  os=1  and 
an=0.1  while  cc  is  varied  over  the  interval 

1 0  <  a  <  1 00  .  (8.43) 

Dependence  of  (SDR)j  on  the  clutter  power  level  is  plotted  in  Fig. 8.23  for  i=I,  II, 
and  m.  Once  again,  case  III  does  not  differ  from  case  n  for  pc=0  due  to  the  flat 
nature  of  the  clutter  power  spectral  density.  The  nonlinearity  is  driven  harder  as  oc 
is  increased.  This  results  in  greater  distortion  of  the  clutter  power  spectral  density 
when  pc=0.9.  Since  the  transversal  filter  for  case  III  is  optimized  to  this  distortion 


(SDR)j  -dB-  (SDR)i  -dB-  (SDR);  -dB- 


149 


(SCR)in 


Figure  8.23  -  Dependence  of  (SDR)m  on  the  Clutter  Power  Level 
m.05<g<0.4975):  ( 1 )  (SDR),  - dB (2)  TSDR)„  -dB-.  (3)  (SORh^-dB-. 
a)  UL=0.5.  p.=0.9.  bffrVf^O.5.  o,=0  cl  fn/L=0~25,  p^0.9, 

(d)  ln/L=0.25,  p^O,  e)  frVfrO.  p>0.9, 7)  o^=0. 

whereas  it  is  designed  on  the  basis  of  a  linear  receiver  for  case  II,  the  output  sig- 
nal-to-disturbance  ratio  is  much  larger  for  case  ID  than  for  case  II  (see  Figs.8.23 

s 

(c)  and  (e)). 

8.3.5  Dependence  of  (SDR)ttt  on  the  Doppler  Frequency 

Plots  of  (SDR)j  versus  fjyfs  are  shown  in  Fig.8.24  for  i=l,  II,  and  in.  The  ca¬ 
ses  of  pc=0.9  and  pc=0  are  considered.  When  pc=0.9,  all  three  curves  of  (SDFty 
have  their  maxima  at  fD=0.5fs.  This  corresponds  to  the  situation  where  the  desired 
signal  spectrum  is  positioned  at  the  minimum  of  the  clutter  power  spectral  density. 
When  pc^O,  the  output  signal-to-disturbance  ratios  are  independent  of  fo/fs  due  to 
the  flat  nature  of  the  disturbance  power  spectral  density.  As  would  be  expected, 
note  that 


(SDR)j>  (SDR)m>  (SDR)n.  (8.44) 
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8.3.6  Dependence  Of  (SDRW  on  the  Number  of  Samples,  N 

Plots  of  (SDR)j,  for  i=I,  II,  and  HI,  are  shown  in  Fig.8.25  for  fc>/fs=0.5,  0.25,  0 
and  pc=0.9,  0.  Once  again,  the  output  signal-to-disturbance  ratios  increase  with 
an  increase  in  the  number  of  samples  processed.  The  amount  of  improvement  is 
approximately  the  same  for  cases  I,  II,  and  in.  As  explained  in  section  8.1 .5,  the 
improvement  is  due  to  1 )  the  power  spectral  density  of  the  desired  signal  becom¬ 
ing  more  concentrated  about  fo  and  2)  the  increase  in  the  degrees  of  freedom  of 
the  transversal  filter  as  N  is  increased. 
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8.4  Comparison  of  Cases  I.  II.  and  III 

Recall  that  (SDR)jn  denotes  the  signal-to-disturbance  ratio  at  the  receiver  in¬ 
put  while  (SDR)j,  where  i=I,  II,  and  III,  denotes  the  signal-to-disturbance  ratio  at 
the  transversal  filter  output.  Since  the  signal-to-disturbance  ratios  are  expressed 
in  dB,  the  processing  gain  for  Case  i,  i=I,  II,  and  m,  is  defined  to  be 

(PG)i=  (SDR)j-  (SDR)jn  .  (8.45) 

As  seen  in  the  previous  sections,  the  nonlinearity  causes  a  reduction  in  (SDR)j  for 
cases  II  and  III  when  compared  to  case  I.  For  i=n,  in,  let  this  loss  be  denoted  by 

(NLL)j  =  (SDR)j-  (SDR) j  .  (8.46) 

Finally,  by  tailoring  the  weights  of  the  transversal  filter  to  account  for  the  nonlin¬ 
earity,  case  III  outperforms  case  II.  The  amount  by  which  (SDR)m  exceeds 
(SDR)n  is  denoted  by 


J  =  (SDR)in-  (SDR)jj 

=  (NLL)n-(NLL)m  .  (8.47) 

Letting  pc=0.9,  N=15,  og=5,  an=0.1,  cs=1.0  and  using  the  numerical  values  given 
in  Eqs.(8.6)-(8.10),  results  were  obtained  for  the  pairs  (a=0.5,  ac=10)  and 
(a=0.05,  ac=100).  These  are  tabulated  in  Tables  8.1  and  8.2,  respectively.  Note 
that  both  Tables  apply  to  the  same  nonlinearity.  However,  the  nonlinearity  for  Ta¬ 
ble  8.2  is  driven  20dB  harder  than  that  of  Table  8.1.  Because  of  the  relatively  nar¬ 
row  power  spectral  density,  the  maximum  processing  gain  is  seen  to  occur  for 
f[)/fs=0.5.  The  loss  in  SDR  caused  by  the  nonlinearity  is  maximum  for  fo/fs=0.5 
and  decreases  as  fQ/fs  approaches  zero.  As  would  be  expected,  (NLL)j  is  much 
larger  for  a=0.05  than  for  a=0.5.  With  respect  to  the  Tables, 

0.39  (dB)  <2.49  (dB)  (8.48) 

for  a=0.5  whereas 
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1.56  (dB)  <  J  £  19.92  (dB)  (8.49) 

for  a=0.05.  The  improvement  of  case  III  over  case  II  tends  to  be  larger  for  smaller 
values  of  the  signal-to-disturbance  ratio  at  the  transversal  filter  output. 
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-dB- 

50.67 

49.96 
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21.46 

9.74 

5.8 
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3.36 

(pG)m 

-dB- 
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28.69 

25.2 

17.87 

8.79 

5.16 

(NLL)|, 

-dB- 

22.07 

21.78 

16.82 

6.08 

3.01 

2.44 

$  #*••••* 

•'  Z.: 
iipi 

a  ; 

-dB- 

21.68 

21.27 

14.51 

3.59 

0.95 

0.64 

J 

-dB- 

0.39 

0.51 

2.31 

2.49 

2.06 

1.8 

Table  8. 1 

Relative  Performance  of  cases  I.  II.  and  III  when  a=0.5. 

In  this  work  we  have  emphasized  the  situation  where  the  signal  is  much  larg¬ 
er  than  the  noise  but  is  much  smaller  than  the  clutter.  Some  general  conclusions 
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-dB- 

26.29 

24.84 

21.59 

15.74 

8.13 

4.75 

(NLL)U 

dB- 

44.27 

41.35 

29.47 

23.55 

21.64 

20.93 

(NLL)UI 

-dB- 

42.71 

37.65 

19.47 

5.94 

1.72 

1.13 

a 

-dB- 

1.56 

3.7 

10.00 

17.61 

19.92 

19.80 

Table  8.2 

Relative  Performance  of  cases  /.  II.  and  III  when  a =0.05. 

that  can  be  drawn  from  the  numerical  results  presented  in  this  Chapter  are  given 
below: 

1)  The  nonlinearity  can  cause  a  significant  reduction  in  the  processing  gain 
of  the  transversal  filter. 

2)  The  signal-to-disturbance  ratios  for  the  three  cases  satisfy  the  inequality 


(SDR)  j  >  (SDR)  IH  £  (SDR)  n  .  (8.50) 
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3)  The  nonlinearity  tends  to  broaden  the  clutter  power  spectral  density.  How¬ 
ever,  a  white  spectrum  at  the  input  of  the  nonlinearity  remains  white  at  its  output. 

4)  (SDR),,  where  i=I,  II,  and  in,  is  maximized  when  the  peak  of  the  target 
signal  spectrum  is  positioned  at  the  minimum  of  the  clutter  spectrum. 

5)  When  the  peak  of  the  target  signal  spectrum  is  offset  from  the  peak  of  the 
clutter  spectrum,  (SDR)i(  where  i=I,  n,  and  in,  tends  to  increase  as  the  clutter  cor¬ 
relation  is  increased  due  to  the  narrowing  of  the  clutter  spectrum. 

6)  The  clutter  model  constrains  the  area  of  the  clutter  power  spectral  density 
to  remain  constant  at  2oc2  as  the  clutter  correlation  is  varied.  Consequently,  the 
clutter  spectral  peak  at  f=0  becomes  larger  as  pc  increases  and  the  spectrum  be¬ 
comes  narrower.  Therefore,  when  the  peak  of  the  target  signal  spectrum  is  cen¬ 
tered  at  f=0,  (SDR),,  where  i=I,  II,  and  III,  tends  to  decrease  with  increasing 
values  of  pc. 

7)  For  pc=0,  the  clutter  is  white  and  the  total  power  spectral  density  of  the 
disturbance  is  a  constant.  Hence,  (SDR)j,  where  i=I,  II,  and  III,  is  independent  of 
the  Doppler  shift,  fD. 

8)  Because  the  shape  of  the  clutter  power  spectral  density  is  modeled  as 
Gaussian,  the  tails  fall  off  exponentially.  Thus,  even  though  the  average  clutter 
power  is  much  larger  than  the  average  noise  power,  the  noise  power  can  limit  per¬ 
formance  when  the  peak  of  the  target  signal  spectrum  is  positioned  in  the  tails  of 
the  clutter  spectrum,  l.  en  an  increase  in  the  noise  power  can  cause  a  significant 
reduction  in  the  output  signal-to-disturbance  ratio.  However,  when  the  target  sig¬ 
nal  spectrum  is  positioned  such  that  the  clutter  limits  performance,  a  moderate  in¬ 
crease  in  the  average  noise  power  has  negligible  effect  on  the  output  signal-to- 
disturbance  ratio. 

9)  When  the  peak  of  the  target  signal  is  positioned  in  the  tails  of  the  clutter 
spectrum  and  the  clutter  spectrum  is  relatively  narrow,  a  modest  increase  in  the 
clutter  power  does  not  affect  (SDR)j,  where  i=I,  n,  and  III.  However,  when  the  tar¬ 
get  signal  spectrum  is  positioned  such  that  the  clutter  is  dominant  over  the  noise, 


157 


(SDR)i  decreases  as  the  clutter  power  is  increased. 

1 0)  As  the  number  of  samples,  N,  grows,  the  output  signal-to-disturbance  ra¬ 
tio  increases  because  (i)  the  power  in  the  target  signal  becomes  more  concentrat¬ 
ed  in  the  frequency  interval  centered  at  fo  and  (ii)  the  degrees  of  freedom  of  the 
transversal  filter  increase. 

11 )  As  a  increases,  the  system  becomes  more  linear  and  the  performance  of 
case  II  approaches  that  of  case  I. 

12)  Because  of  the  widening  by  the  nonlinearity  of  the  clutter  power  spectral 
density,  the  output  signal-to-disturbance  ratio  becomes  less  sensitive  to  an  in¬ 
crease  in  the  average  noise  power.  However,  as  explained  in  section  8.2.3,  the 
output  signal-to  disturbance  ratio  can  increase  with  an  increase  in  noise  power. 

13)  When  the  clutter  is  white,  the  decrease  in  output  signal-to-disturbance 
ratio  is  identical  for  cases  I  and  II  as  the  average  clutter  power  is  increased.  How¬ 
ever,  when  the  clutter  is  not  white,  the  degradation  is  greater  for  case  n  than  for 
case  I. 

14)  Case  III  is  better  able  to  outperform  case  II  when  the  clutter  is  more  cor¬ 
related  (i.e.,  the  clutter  spectrum  is  narrower)  because  the  nonlinearity  causes  a 
greater  distortion  in  the  clutter  spectrum.  Cases  II  and  in  perform  the  same  when 
the  clutter  is  white. 

1 5)  When  the  clutter  is  not  white,  case  III  performs  increasingly  better  rela¬ 
tive  to  case  n  as  the  nonlinearity  is  driven  harder. 

16)  For  a  narrow  clutter  power  spectral  density,  case  in  performs  increasing¬ 
ly  better  relative  to  case  II  as  the  average  noise  power  is  reduced. 

17)  As  can  be  seen  from  Fig. 4. 4,  the  output  of  the  nonlinearity  is  decidedly 
nonGaussian  when  a<1 .  Yet,  the  Gaussian  receiver  is  surprisingly  robust  with  re¬ 
spect  to  the  signal-to-disturbance  ratio  at  the  output  of  the  transversal  filter,  espe¬ 
cially  when  a=l. 

18)  When  the  receiver  is  nonlinear,  use  of  the  transversal  filter  optimized  for 
the  linear  receiver  can  result  in  (SDR)n  being  smaller  than  (SDR)jn,  as  shown  in 
Table  8.2. 
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Appendix  A 

Derivation  of  Equation  (4.3) 

In  general,  the  received  signal  r^t)  will  be  a  modulated  signal  with  carrier 
frequency  fc.  The  quadrature  components  of  r^t)  are  given  by 


r,c(t)  =  [r,<t).cos(2*fct)]|Lp 

(A- 1 ) 

^  r,s(,)  =  [ri(,)sin(2lt,ct)1,Lp 

(A.2) 

where  [.]|Lp  denotes  the  operation  of  passing  the  argument  through  an  ideal  low 
pass  filter.  The  received  signal  can  then  be  expressed  as 

ri(t)  =  r1c(t)cos(27tfct) -rls(t)sin(2itfct) 

j2ftfCt  ,a  OX 

=  Re{[r1c(t)+jrls(t)]e  c}  •  <A3) 

However,  from  Eq.(3.4),  r^t)  is  related  to  its  complex  envelope  according  to  the 
relation 

j27tf  t 

r1  (t)  =  Re  {?,  (t)  •  e  c  }  -  (A.4) 


Comparison  of  Eqs.(A.3)  and  (A.4)  reveals  that 

Mt)  *  riC(t)  +jris(t) 

Use  of  Eqs.(A.I)  and  (A.2)  in  Eq.(A.5)  results  in 

j2*  U 

rt  (t)  -  {rt  (t)  -  e  c}  * 

1  ILP 


(A.5) 


(A.6) 


This  operation  is  depicted  in  Fig.A.1  where  gn_p(t)  denotes  the  impulse  response 
of  the  ideal  low  pass  filter  whose  transfer  function  is  shown  in  Fig.A.2.  The  band 
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width  W  is  chosen  equal  to  half  the  bandwidth  of  r^t). 

From  Fig.A.1,  it  is  obvious  that  the  complex  envelope  is  determined  by  the 
convolution  of  r1(t)d2n,ct  with  gn_p(t).  It  follows  that 

7  i2nfcx 

Mt)  =  J  r1  (x)  •  e  c  9,Lp(t-x)dx  .  (A.7) 

— o© 
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As  a  result 


j2rcfcx, 


E[*i  (t1)r1  (t2) ]  =  E  J  r,  (x.,)  •  e  c  ■  -x1)dx1 


I  r1  <x2> 


j2refcx2 


^ILP  ^2  ~  x2^  dx2 


*  j  J  E[^  (x1 )  ^(x^Je 


j27tf  (X1+X2) 


*  9|lp  (^  ~  xi)  ‘  9(|_p  (^2  ~  x2^xi^x2  ■  (A.8) 

Assuming  the  received  waveform  from  a  given  range  cell  is  stationary,  the  auto¬ 
correlation  function  has  the  property  that 


R^x^xg)  =  E  [r1  (x.,)  r1  (x2) ]  =  R1(x1-x2) 


(A.9) 


Hence,  Eq.(A.8)  becomes 

oo  oo 

E  [r-j  (t1 )  rn  (t2)  ]  =  \  1R.  (xi~x2)e 


j27lfc(X1+X2) 


9|LP(ti  xi^ 


9|LP(t2-X2)dxldx2 


(A.  10) 


In  terms  of  the  power  spectral  density  of  r-j(t) ,  which  is  denoted  by  Si(f),  the  auto¬ 
correlation  function  can  be  expressed  as 


R1(x1-x2)=  JSl(f).ej27lf(Xl  Xs)df  . 


(A.  11) 


Substitution  of  Eq.(A.II)  into  Eq.(A.IO)  yields 

EIT, (.,)?,(.,)]  =  J  J  Js,(0  ei^x.-x*>  + 


•  9|LP<'l-X1)'9||p02-><2>d><ldx2d< 
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oo  oo 

=  J  (f)  ■  J  9j|_p (t-i  ~xi) 

— OO  — oo 


j2KX1(f  +  f  ) 
e  c  dx1 


7  -j2*X2(f-f  ) 

"  J  9||_P  ^2  ~  X2^  ‘  ®  ^X2 


7  [7  -j2*vl(uf) 

J  si(f)  j  J  g,LP(vi)  e  dv1 

— OO  V-OO 


oo 

J 


gii_p(v2) 


j27lV2(f-f  ) 

e  c  dv2 


j2Jlt1(f  +  f  )-j2Jlt2(f-f  ) 
e  c  c  df 


(A.  12) 


However, 


G|Lp  (0  ~  J  9|(_p  (1)  *®  dt 


(A.  13) 


Use  of  Eq.(A.13)  in  Eq.(A.12)  yields 

oo 

E  [r-j  (t-| )  r^  (t2)  ]  =  J  S.|  (f)  •  G|j_p  (f  +  f^,)  ■  ^ilp  ~ fg) 


j27lt1(f  +  f  )-j27lt2(f-f_) 
e  c  c  df  .  (A.  14) 


Assuming  the  modulation  of  r-|(t)  to  vary  slowly  relative  to  variations  of  the  carrier, 
fc  »  W  .  It  follows  that  the  transfer  functions  GtLp(f+fc)  and  GiLP*(f~fc)  wil1  not 
overlap  in  frequency.  Thus, 

GILP(f +  fc)  GILP*  (f"fc>  =  0  '  (A.  15) 


Consequently,  the  integrand  in  Eq.(A.14)  is  zero  and 
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E  [r,  (t1)?1  (t2)]  =  0  .  (A.  16) 

When  ti  » t  and  tz  =  t-x,  this  result  becomes 

E  [r.,  (t)  r1  (t  —  x)  ]  =  0  .  (A.  17) 
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Appendix  B 

Derivation  of  Equations  (4.56)  and  (4.57) 


Recall  from  Eq.(3.30)  that 


2L 

©o 

1+ 1 
k.  i 

f-l'l 

k  n2k  +  i 

ck  (2k-  1)!!  f  2k 

fc  kite2 

V  g 

o2 

^  QJ 

22k  (2k  -Hi)!  I1! 

ri  .  (B.1) 


To  simplify  this  expression,  let 

G  = 


2L 


■ 

t 

\ 


2no£ 


and 


'-i' 


k  r*  2k  +  1 


(2k-  1 ) !! 


^02j  22k  (2k  +  1 ) ! 


(B.2) 


(B.3) 


Eq.(B.1 )  then  becomes 


=  G 


+ 


I* 

k  =  i 


-  ,2k  - 
k'lrtl  ‘ri 


(B-4) 


As  seen  in  Eq.(4.55),  the  relation  for  R3(t,x)  involves  E 


r2(t)  f  r2(t-x) 


cJ 


and 


r2(t) 


* 


r2  (t-x) 


cJ 


.  The  evaluation  of  these  two  expectations  is  considered 


in  this  appendix. 

We  first  consider 


r2  (t)  r2(t-x) 

rc 


f 


=  G2E 


2k 

?i  (t)  +  x  AkTri  (t)l  "ri  (<) 

k  =  i 


?1  (t-x)  +  X  A  -(?!  (t-x)|2p- ?!  (t-X) 

P=1 
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*  G2E  r,  (t) rt  (t-x)  +  £  Ap  rt  (t)  •  [r1  (t-t)]2p  +  1  •  [?/  (t-x) J2p 
L  P=i 

oo 

+  I  Ak-  [r,  (t))2k-1  -  [?,'  (t)l2k-r,  (t-x) 
k  =  i 

+  X  £  AkA  •rr1(t)]2kt,  .[r1*(t)]2k.(-r1(t-T)l2P+1 

k  =  ip  =  i  p 

•  [r,'(t-x)l2p]  .  (B.5) 

Interchanging  the  order  of  summation  and  expectation,  we  obtain 
E  r2 (t)  ?2 (t-x)  =  G2E [r1  (t)  r,  (t  -  x)  ] 

c  TcJ 

+  G2  X  ApE^1  (t-x)]2p  +  1  •  [?/  (t  —  X)  ] 2p) 

P  =  1 

+  G2  I  AkE[[?1(t)]2k  +  1-  [?/  (t)]2k  F1  (t-x)" 

k  =  i 

+  q2S  I  AkA  E[[?l(t)]2k+1.  [?;  (t)]2k.  [?l(t-x)]2p+1 

k= ip  =  i  p 

•  lK  (t-x)]2p]  .  (B.6) 

From  Eq.(A.17) 


E  [r,  (t)  ?!  (t-x)l  =  0  . 


(B.7) 


Hence,  the  first  term  on  the  right  hand  side  of  Eq.(B.6)  is  zero. 

To  proceed  further,  recall  that  r1  (t)  is  modeled  as  a  zero  mean  stationary 

complex  Gaussian  random  process.  If  z-|,z2 . z2rn  denote  2m  jointly  Gaussian 

random  variables,  it  is  known  [7]  that 

E[2122 . 22m>  ■  I 

all 

pairs 

Each  term  in  the  summation  of  Eq.(B.8)  contains  the  product  of  m  pair  wise  ex¬ 
pectations  where  each  of  the  2m  random  variables  appears  once  and  only  once. 
The  number  of  terms  in  the  summation  equals  the  number  of  different  ways  2m 

different  variables  z-|,z2 . z2m  can  be  chosen  in  pairs,  which  is 

(2m) !/  (2mm!) .  For  example,  when  m  =  2,  the  summation  in  Eq.(B.8)  consists 
of  three  terms  as  given  by 

E  [z1z2z3z4]  =  E  [z,z2]  E  [z3z4]  +  E  [z1z3]  E  [z2z4] 

+  E[z1z4]E[z2z3]  .  (B.9) 

Now  examine  E  (r1  (t)  [r1  (t-t)]2p  +  1  [r1  (t-x)]2P}.  This  expectation 
involves  the  product  of  (4p+2)  jointly  Gaussian  random  variables.  Consequently,  a 
typical  term  in  Eq.(B.8)  consists  of  a  product  of  (2p+1)  pair  wise  expectations.  To 
avoid  a  zero  result,  it  is  necessary  that  at  least  one  of  these  terms  not  have 
E  [rt  (t)  rt  (t  —  x)  ]  as  a  factor.  However,  this  is  not  possible.  The  term  with  the 
fewest  number  of  factors  of  the  form  E  [r1  (t)  -  r1  (t  —  x)  ]  arises  when  r1  (t  -  x)  is 
paired  with  r^  (t-x)  a  total  of  2p  times  and  r1  (t)  is  paired  with  r1  (t-x)  once. 
All  other  terms  in  the  summation  of  Eq.(B.8)  have  factors  in  which  r1  (t)  is  paired 
with  r1  (t  -  x)  more  than  once.  It  follows  that 
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Substitution  of  these  results  into  Eq.(B.6)  yields  the  conclusion  that 


E  r2(t)  f  r2  (t  -  x)  f  =  0  . 
c  'c. 


(B-14) 


We  next  consider  E 


r2  (t)  r2  (t  -  x)  .  From  Eq.(B.4)  it  follows  that 
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G2E  [r1  (t)V  (t  —  x)  ] 


.  i. 


+  Qi£  I  ApEt'r^t)  •  [rt  (t-x)]2p  +  1-  [r1(t-x)]2p} 

p  =  i 

+  G2  £  AkE[[r1(t)]2k+1.  (V  (t)]2k  f/  (t-x)] 
k  =  i 


oo  oo 


2k  *  ..  .  _2P  +  1 


+  G2  £  £  AkA  E[(f,(t))2k  +  1-  [?,*  (t)l  -  [r,  (t-x)] 

k»  ip  =  i  p 


•  [?1<t-x)]2p]  . 


(B.15) 


Each  expectation  can  be  expanded  as  a  sum  of  products  of  pair  wise  expecta¬ 
tions  as  given  by  Eq.(B.8).  In  each  summation  it  is  readily  seen  that  there  exists  at 
least  one  term  that  is  nonzero.  The  nonzero  terms  involves  factors  of  the  form 


-  * 


R,(0)  =  E[r,(t)r,  (t)  J 


R,  (xj  =  E [r,  (Of,  (t-x)] 


_  * 


* 


Ri  (0)  =  E^  (t-x) r.,  (t-x)] 


(B.16) 


«  * 


* 


Rt  (x)  =  E  [r,  (t)  r,  (t - x)  ]  . 


From  this  and  assuming  stationarity  of  (t) ,  it  is  readily  concluded  that  the  auto¬ 


correlation  E 


r2(t) 


r2  (t-x) 


'cj 


is  a  function  of  x  alone. 
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